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Multi-objective optimization of high speed 
vehicle-passenger catamaran by genetic 
algorithm 


Part ll 
Analysis of the results 


Zbigniew Sekulski, Ph. D. 
West Pomeranian University of Technology in Szczecin 


ABSTRACT 


Real ship structural design problems are usually characterized by presence of many 
conflicting objectives. Simultaneously, a complete definition of the optimum structural 
design requires a formulation of size-topology-shape-material optimization task unifying 
the optimization problems from the four areas and giving an effective solution of the 
problem. Any significant progress towards solving the problem has not been obtained so 
far. An objective of the present paper was to develop an evolutionary algorithm for multi- 
objective optimization of the structural elements of large spatial sections of ships. Selected 


elements of the multi-criteria optimization theory have been presented in detail. Methods for solution of the 
multi-criteria optimization problems have been discussed with the focus on the evolutionary optimization 
algorithms. In the paper an evolutionary algorithm where selection takes place based on the aggregated 
objective function combined with domination attributes as well as distance to the asymptotic solution, is 
proposed and applied to solve the problem of optimizing structural elements with respect to their weight 
and surface area for a high - speed vehicle-passenger catamaran structure, with taking into account several 
design variables such as plate thickness, scantlings of longitudinal stiffeners and transverse frames, and 
spacing between longitudinal and transversal members. Details of the computational models were kept at 
the level typical for conceptual design stage. Scantlings were analyzed by using the selected classification 
society rules. The results of numerical experiments with the use of the developed algorithm are presented. 
They show that the proposed genetic algorithm may be considered an efficient tool for multi-objective 
optimization of ship structures. 
The paper has been published in the three parts: Part I: Theoretical background on evolutionary multi- 
objective optimization, Part II: Computational simulations, and Part III: Analysis of the results. 


Keywords: ship structure; multi-objective optimization; evolutionary algorithm; 
genetic algorithm; Pareto domination, set of non-dominated solutions 


ANALYSIS OF THE RESULTS AND 
CONCLUSIONS DRAWN FROM THE 
COMPUTATIONAL SIMULATIONS 


Three series of the computer simulations, signed sym1, 
sym2 and sym3, confirmed effectiveness of the developed 
computational algorithm and computer code for solution 
of the formulated problem of ship structure topology-size 
multi-objective optimization. As a result of the calculations 
an approximation of the Pareto-optimum set containing, in 
each simulation, from a few to more than ten non-dominated 
solutions, was found. The obtained results do not allow to 
unequivocally conclude which of the examined factors: (1) 
objective function aggregation strategies, (2) domination 
attributes included into selection process, and (3) distance to 


asymptotic solution included into selection process, is most 
advantageous. 

In the case of studying the influence of optimization 
criteria aggregation strategy, visual assessment of the shape 
of the obtained approximations of the Pareto-optimum set 
suggests an advantage of the strategy with random values of 
the weight coefficients w, (sym1-2) and the least effectiveness 
of the strategy with fixed values of the weight coefficients w, 
(sym1-1). 

The effectiveness of the strategy with random selection of 
single optimization criteria in the selection process (sym1-3) 
is intermediate. In the case of the constrained problems it also 
turns out that the components of the penalty functions introduce 
a random contribution to the fitness function thus causing the 
strategy with the fixed weight coefficients w, to be practically 
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a strategy similar to the two others and it also allows to find 
approximation of the Pareto-optimum set with the adequate 
accuracy. 

From the found sets of the compromise solutions a user can 
select, in the next stage, one or a few solutions by applying 
additional premises which are not included in the optimization 
model. He can also select suggested non-dominated solutions 
the closest to the asymptotic solutions f=: 


f [1086.28 7422.10]™[kKN m2] => 


*syml-1 = 


= fi symı-1X) = 1086.28 KN, Ê symı-1X) = 7422.10 m? 


P ymt-2 = [1113.65 7361.45]™ [kN m?] = 
=> fh sym-2(X) = 1113.65 KN, f, symı-2(X) = 7361.45 m? 
Fmi-3 = [1153.68 7381.57]" [KN m?] = 


= fi gymt-a(X) = 1153.68 KN, fy sm1- 3(¥) = 7381.57 m? 


a) 7800 
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Since this way three solutions are obtained, the next 
question is which of them can be recommended as the best». 
Here the following procedure is suggested by this author: non- 
dominated solution sets obtained in subsequent simulations can 
be merged into a temporary solution set presented in Fig. 39a. 
In this set only a part of solutions is non-dominated ones, 
Fig. 39b. In the set of 15 non-dominated solutions obtained 
by using the results of three simulations, a distance of each of 
them from the asymptotic objective in normalized objective 
space, can be determined, Fig. 39c. The least distance equal 
to 1.082, was obtained for the solution f,(x) = 1113.65 kN and 
f,(x) = 7361.45 m? found in the simulation sym1-2 (random 
values of the weight coefficients w, and w, in the range [0, 1 ]). 
The solution can be recommended as a single solution of the 
formulated problem of multi-objective optimization. 

Effectiveness of the three examined multi-objective 
optimization strategies which use optimization criteria values 
and functions representing violation degree of constraints, 
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Fig. 39. Selection of single, recommended solution of multi-objective optimization problem by using non-dominated solution sets obtained in the three series 
of computer simulations: sym1-1, sym1-2, sym1-3: a) temporary set composed of non-dominated solutions of each simulation, b) selection of non-dominated 
solutions in temporary set, c) determination of distance of non-dominated solutions from asymptotic objective in normalized objective space, 

d) values of optimization criteria for the found closest solution: f,(x) = 1113.65 kN and f,(x) = 7361.45 m 


D Let us remember that in the multi-objective optimization there is not the single best solution of the problem and the formulated 
recommendation should be treated as a subjective choice by a person taking decision. 
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can be roughly evaluated also by assuming that the number 
of solutions produced by a given algorithm, that is a part 
of set of non-dominated solutions obtained on the basis of 
results produced by all algorithms, Fig. 39, is a measure of 
algorithm effectiveness. In the examined example the particular 
simulations: syml-1, syml-2 and sym1-3 have brought 
respectively 6, 6 and 3 non-dominated solutions into set of 
non-dominated solutions. On this basis a conclusion may be 
suggested that the strategies of fixed values of the optimization 
criteria weight coefficients w, = w, = 0.5 and random values 
of the weight coefficients w, and w, in the range [0, 1] show 
similar effectiveness whereas the strategy of random values of 
the weight coefficients w, and w, equal to 0 or 1 proved to be 
the least effective. The issue of examining the effectiveness of 
multi-objective evolutionary algorithm is very relevant and not 
fully solved hence it requires a separate research which exceeds 
however content of this case study. 

The conducted analysis of computer simulation results 
of the problem of ship structure multi-objective optimization 
in question allows to state that in the studied cases the most 
effective strategies were the following: (1) that with random 
values of the weight coefficients w, and w, in the range [0, 1], 
and that with fixed values of the optimization criteria weight 
coefficients w; = w, = 0.5. Less effective was the strategy with 
random values of the weight coefficients w, and w, equal to 
0 orl. 

The recommended non-dominated solution was obtained 
for the values of design variables specified in Tab. 8. The 
corresponding dimensions of the ship cross-section are given 
in Fig. 40. 

From the study of influence of dominance attributes and 
distance from asymptotic solution on the effectiveness of the 
algorithms of sym2-1, sym2-2 and sym2-3, also satisfactory 
results were achieved as follows: 


P m21 = (1105.95 7345.11]™ [kN m2] = 
> fi gymo-1(X) = 1105.95 KN, f mai) = 7345.11 m? 


Pym- = [1192.04 7327.41 ]™ [kN m?] = 


=> f gymna-2(X) = 1192.04 KN, fy yma) = 7327.41 m? 


f m23 = [1060.03 7485.93]"-[kN m2] > 
= fy „m2.3(X) = 1060.03 KN, fyyna-3(X) = 7485.93 me 


In the case of necessity to identify a single solution from 
a series of simulations one can apply the earlier described 
procedure of aggregation of set of non-dominated solutions 
obtained in particular simulations and identify the non- 
dominated solution nearest the asymptotic solution. 

The performed calculation investigations have positively 
verified the effectiveness of the combined fitness multi- 
objective evolutionary algorithm developed by this author, as 
well as the calculation tool built for solving the unified ship 
structure topology-size multi-objective optimization problem. 
Particular computer simulations have produced a dozen or 
somewhat more of non-dominated solutions which constitute 
the set of trade-off solutions from among which decision makers 
may choose one or more of them for further development. The 
algorithm developed as a part of the underlying work allows 
also to pinpoint a single variant closest to the asymptotic 
solution which may be proposed as a single solution of the 
multi-objective optimization problem. 

Fig. 41 presents the comprehensive results of the multi- 
objective optimization of the ship hull structure: (1) general 
arrangement and main particulars of an example ship, (2) 


optimization criteria, (3) simulation main parameters and 
control variables, (4) values of optimization criteria for the 
obtained non-dominated variants, (5) values of optimization 
criteria for the variant closest to the asymptotic solution, and 
(6) structural dimensions and scantlings for this variant. 


Tab. 8. Values of design variables recommended 
as a result of multi-objective optimization 


Description 


serial No. of mezzanine deck plate 


serial No. of mezzanine deck bulb 


serial No. of mezzanine deck T-bulb 


number of web frames 


number of mezzanine deck stiffeners 


serial No. of superstructure I plate 


serial No. of superstructure I bulb 


serial No. of superstructure I T-bulb 


SOS) Ola ni ny; RAJU Nj] = 


number of superstructure I stiffeners 


serial No. of inner side plate 


serial No. of inner side bulb 


serial No. of inner side T-bulb 


number of inner side stiffeners 


serial No. of bottom plate 


serial No. of bottom bulb 
serial No. of bottom T-bulb 


number of bottom stiffeners 


serial No. of outer side plate 


serial No. of outer side bulb 


serial No. of outer side T-bulb 


number of outer side stiffeners 


serial No. of wet deck plate 
serial No. of wet deck bulb 
serial No. of wet deck T-bulb 


number of wet deck stiffeners 


serial No. of main deck plate 
serial No. of main deck bulb 
serial No. of main deck T-bulb 


number of main deck stiffeners 


serial No. of superstructure II plate 


serial No. of superstructure II bulb 


serial No. of superstructure II T-bulb 


number of superstructure II stiffeners 


serial No. of upper deck plate 


serial No. of upper deck bulb 


serial No. of upper deck T-bulb 


number of upper deck stiffeners 
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Frame spacing: 1250 mm 


Fig. 40. Ship structural dimensions and scantlings recommended as a result of multi-objective optimization; 
structural material: for plates - 5083-H111 Al alloy, for profiles - 6082-T6 Al alloy 


PERFORMANCE ASSESSMENT OF 
THE DEVELOPED MULTI-OBJECTIVE 
OPTIMIZATION ALGORITHM 


The question should be answered if the proposed combined 
fitness multi-objective evolutionary algorithm (Eq. 28) is more 
efficient than the evolutionary algorithms used in the process 
of selection of only optimization criteria (in the form of scalar 
substitute optimization criteria) and functions representing the 
degree of constraint violation (Eq. 8). Unfortunately the answer 
cannot be simple and unequivocal. 

In practical problems the ship structural multi-objective 
optimization which produces a set of Pareto-optimum 
solutions may be very computation time-consuming or even 
impossible to be performed. In the cases the evolutionary 
algorithms of multi-objective optimization do not guarantee 
identification of Pareto-optimum compromises but can 
help identifying a satisfactory approximation, i.e. a set of 
solutions hoped to be not too far distant from searched front 
of optimum solutions (in the sense of Pareto). However in this 
case methods are necessary to evaluate how good produced 
solutions of formulated problems are. And, this leads to 
the question: how to compare effectiveness of different 
algorithms? In the context of the presented work the question 
may be formulated as follows: how to compare effectiveness 
of the studied evolutionary multi-objective optimization of 
ship hull structure, assumed for different evaluation strategies 
of fitness function. 

In the case of the evolutionary algorithms of multi- 
objective optimization, statistical in their nature, evaluation 
of obtained results and comparison of effectiveness of 
optimization algorithms implementing different strategies 
is a very difficult task, arousing much controversy and 
misunderstanding. Whereas visual and qualitative comparison 
of the sets approximating Pareto front is commonly used 
for deduction of quality of evolutionary multi-objective 
optimization, in the case of quantitative methods the searching 
of proper standards are just under way [Knowles, Thiele, 
Zitzler (2006)]. 

Generally accepted procedures enabling to compare quality 
of solutions obtained by different algorithms (usually in many 
runs) or solutions of the same algorithm obtained in many runs, 
quantitatively taking into account their statistical characteristics 
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are necessary [Sarker, Coello Coello (2002)], [Knowles, Thiele, 
Zitzler (2006)]. 

The problem is extremely difficult also because in this case, 
as opposed to single-objective optimization, it is necessary to 
compare not individual solutions, but vectors representing sets 
of non-dominated solutions under assumption that they are an 
approximation ofa practically unknown set of Pareto-optimum 
solutions, referred to as a front of Pareto-optimum solutions. In 
consequence, most comparative studies are based on different 
methodologies and assumptions, and therefore results obtained 
from such studies are difficult to be used in mutual comparisons 
[Knowles, Thiele, Zitzler (2006)]. 

In the case of single-objective optimization the selection of 
quality measure is obvious and simple: optimization criterion. 
The quantity is unequivocally defined, optimization criterion 
value calculated for every test solution, and smaller or greater 
depending on the task, which corresponds to better solution. 
In the case of multi-objective optimization it is not clear what 
better” means: is it that located closer to the front of optimum 
solutions, covering a wider range of solution characteristics, 
or something else? And, it should be realized that the front 
of Pareto-optimum solutions is unknown. That is why it is 
difficult to define proper quality measure for approximations of 
an unknown front of Pareto-optimum solutions. Therefore for 
comparing and evaluating results of qualitative evaluation of 
multi-objective evolutionary algorithms, graphical presentation 
of obtained non-dominated solutions has been first of all used 
until recently [Veldhuizen (1999)]. 

In recent years in this field a certain progress has been made 
and several papers concerning the quantitative comparing of 
different approximations of Pareto-optimum set can be found. 
The most popular is the unary quality measure, i.e. that which 
attributes, to every single approximation set, one numerical 
value which reflects a specified quality aspect [Veldhuizen, 
Lamont (2000)], [Zitzler, Thiele, Laumanns, Fonseca, Grunert 
da Fonseca (2002)]. To increase the deducing strength the unary 
quality measures are usually used jointly, to cope of taking 
into account different aspects of the notion of “quality”. Other 
methods are based on binary quality measures which assign 
numerical values to pairs of solutions, [Zitzler, Thiele (1998)], 
[Hansen, Jaszkiewicz (1998)]. 

Third group of evaluation methods, completely different 
in conceptual respect, is the method of attainment function, 


High-speed vehicle-passenger catamaran, 
Auto Express 82 type 


Main particulars: 

- length overall, L = 90.0 m, 

- breadth, B = 23.0 m, 

- maximum depth, H = 11.7 m, 

- length of midship block-section 
subject to optimisation, 17.0 m. 


ay 


nT 


me 


Properties of structural material: 

- aluminium alloys, 

- 5083-H111 alloy for plates (Ro.2 = 125 N/mm’), 

- 6082-T6 alloy for bulb extrusions (Ro2 = 250 N/mm’). 


Multi-objective topology and size optimisation of ship structure 
by evolutionary algorithm with combined fitness function, CFMOEA, 
based on genetic algorithm 


Combined fitness function: 


A(X) = wa f(x) + We fo(X) + Wrank FANK(X) + Weount COUNL(X) + Weistance distance(x) + Y (w - penalty,.(x)) 


k=l 
Optimization criteria: Non-dominated solutions: 
- structural weight f,(x), KN, (1) fy(x1) = 1069 KN, f(x) = 7790 m?, 
- surface area for cleaning and painting f(x), m’. (2) A(x2) = 1077 KN, f(Xx2) = 7742 m?, 
(3) A(x3) = 1085 KN, f(x3) = 7518 mê, 
Parameters of CFMOEA simulation: (4) (x4) = 1106 KN, fo(x4) = 7345 m?, 
- sym2-1, (5) A(xs) = 1316 KN, (xs) = 7231 m?, 
- number of individuals n; = 5,000, (6) A(xs) = 1391 KN, f(xs) = 7130 m°, 
- number of generations ng = 10,000, (7) f)(x7) = 1607 KN, f2(x7) = 7074 m°, 
- W4 = W2 = 0.0, (8) f,(Xg) = 1767 KN, f(x) = 7074 m°, 
- Wrank = 3.0, Weount = 0.0, Waistance = 0.0. (9) (Xo) = 1968 KN, f(x) = 7044 m°, 


(10) f,(x10) = 2244 KN, f(x1o) = 7022 m°. 
The best (recommended) solution: 
- structural weight f,(x4) = 1106 KN, 
- surface area f(x4) = 7345 m’. 


Midship section structural arrangements and scantlings of the best solution 


Frame spacing: 679 mm. 
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Fig. 41. Result of the multi-objective evolutionary optimization of the ship structure with respect to the structural weight f, 
and surface area of structural members for maintenance (cleaning and painting) f, (svm2-1) 
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Fig. 41 cont. Result of the multi-objective evolutionary optimization of the ship structure with respect to the structural weight f, and surface area 


of structural members for maintenance (cleaning and painting) f, (sym2-1), which contains: evolution of the highest value of fitness function 
and the lowest value of non-dominated solution distance from asymptotic one; structure of the set of non-dominated solutions; number 
of generations in which particular non-dominated solutions were found; detailed structure of the set of non-dominated solutions; 
structure of the set of non-dominated solutions in normalized objective space 
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[Fonseca, Fleming (1996)], in which the probability of 
achieving any chosen objectives in evaluation space is assessed 
on the basis of the knowledge of many approximation sets. 

Besides the above mentioned diversity of methods it is still 
unclear what mutual relations of certain quality measures are 
like (what is their mutual connection) and what their relative 
advantages and disadvantages are [Zitzler, Thiele, Laumanns, 
Fonseca, Grunert da Fonseca (2002)]. In consequence there are 
no agreed opinions stating which quality measure or measures 
should be used in specific cases [Zitzler, Thiele, Laumanns, 
Fonseca, Grunert da Fonseca (2002)]. 

To define a reliable evaluation methodology is very 
important for its application to algorithm validation. However, 
as far as the problem of multi-objective optimization is 
concerned there are several reasons due to which it is difficult 
to evaluate obtained results. Firstly, from evolutionary 
algorithms many solutions are obtained instead of only one, 
usually as many as possible solutions belonging to the set of 
non-dominated solutions, approximation set of Pareto-optimum 
solutions, are aimed at. Secondly, evolutionary algorithms 
are stochastic, therefore for effectiveness evaluation it is 
necessary to run many simulations and subject obtained results 
to statistical analysis; in this case drawn conclusions will also 
have stochastic characteristics. Thirdly, we may be interested 
in measuring different aspects of quality; for example, we 
may be more interested in possessing a robust, but slower, 
algorithm convergent to Pareto front practically in every case, 
than in a faster algorithm but convergent to Pareto front only 
occasionally (in case of specific types of tasks); we also may 
be interested in evaluating behaviour of evolutionary algorithm 
in the course of simulation, trying to determine its capability 
of maintaining diversity and gradual convergence to set of 
solutions close to Pareto front. This short discussion shows 
how difficult is to develop effectiveness measures for multi- 
objective optimization evolutionary algorithms. 

The next problem in the discussion is the question: what 
should be measured? It is very important to determine what 
kind of results will be subjected to measurement, evaluation and 
analysis, and to define quality measures according to a task. 

It is obvious that in the formulating of a good quality 
measure for multi-objective optimization evolutionary 
algorithm the following should be considered [Zitzler, Deb, 
Thiele (1999)]: 

1. the minimizing of distance between approximation of Pareto 
front, obtained by the algorithm, and a real Pareto-optimum 
front of solutions, (of course if so is known, what in practice 
of optimizing engineering objects does not happen), 

2. the maximizing of diversity of obtained non-dominated 
solutions, i.e. the arranging of non-dominated solutions 
in approximation set over empirical compromise area, as 
smoothly and homogeneously as possible, 

3. the maximizing of number of solutions in the set which 
approximates the Pareto-optimum set. [Zitzler, Thiele, 
Laumanns, Fonseca, Grunert da Fonseca (2002)], [Zitzler, 
Laumanns, Bleuler (2002)], [Zitzler, Thiele, Laumanns, 
Fonseca, Grunert da Fonseca (2003)], [Fonseca, Knowles, 
Thiele, Zitzler (2005)], [Knowles, Thiele, Zitzler (2006)] 
presented the most extensive (in the author’s opinion) 
review of problems related to evaluation of effectiveness 
of randomized multi-objective optimization algorithms. 
Assuming that a set of incomparable solutions (called 
an approximation set), is a result of operation of multi- 
objective optimization evolutionary algorithm, they 


proposed a mathematical basis for studying multi-objective 
optimization effectiveness algorithms. 


In particular [Zitzler, Thiele, Laumanns, Fonseca, Grunert 
da Fonseca (2002)] showed that if we have two sets of solutions, 
a and B, which approximate Pareto-optimum set of solutions, 
then we cannot elaborate a finite set of quality measures, which 
can indicate if the approximation a of Pareto-optimum set is 
better that the approximation B, in every case. Elaborated 
quality measures may be applied only to specific aspects of 
the quality concept, therefore the only thing left to deduct is 
that the approximation a is not worse that the approximation 
B, which means that either the approximation a is better 
than the approximation B, or the approximations a and B are 
incomparable with regard to a specified quality measure. This 
statement cannot be generalized in a way indicating that the 
approximation a is always better than the approximation B. 
So if it is impossible to state in a close and exact quantitative 
way the supremacy of one approximation set over the other, 
therefore it is impossible to state the supremacy of one 
algorithm over the other. Choice of one of them is determined 
by the efficiency in every specific case it was used. 

[Sarker, Coello Coello (2002)] made a review of propositions 
they considered the most important and enabling to measure 
the three above listed aspects of the notion ,,quality”, subjected 
to evaluation. They also observed that there is no method that 
would allow to measure the three aspects with one value only. 
Unfortunately, attempts to elaborate a single measure for 
grasping them together have not been successful so far because 
they concern very different algorithm characteristics. That is 
why the attempts to reduce them to one measure may lead to 
misunderstandings. Therefore the using of different quality 
measures to estimate different aspects of algorithm behaviour 
seems more proper in practice). 

On the basis of visual assessment of approximating sets 
attained in particular simulations, Fig. 42, it can be stated 
that the best solutions, i.e. approximation sets, were obtained 
in the simulations: sym2-1 (w, = w, = 0.0, Wax = 3.0, 
Weount = 9-0, Woistance = 9-0), sym1-2 (w,, w, are random in the 
range [0, 1], Wan = 9-0, Weount = 0-0, Waistance = 0-0) and sym2-3 
(W, = W2 = 0.0, Wank = 9-0, Weount = 9-0, Waistance = 3-0). Solutions 
obtained from the simulations are arranged the most uniformly, 
hence they can be expected to be a good representation of 
searched Pareto-optimum front. In the case of the remaining 
simulations: syml-1 (w, = W = 0.5, Wax = 9-0, Weount = 0.0, 
Wdistance = 9-0), sym1-3 (w,, w, are random 0 or 1, Wank = 0.0, 
Weount = 9-0, Woistance = 9-0) and sym2-2 (w, = w, = 0.0, Wan = 9.0, 
Weount = 3-9, Woistance = 9-0), the found solutions are arranged 
less uniformly, hence they represent Pareto-optimum front 
in a worse manner. Therefore it can be approximately stated 
that out of six conducted research simulations the following 
were found more effective: the simulation (1) that took into 
consideration influence of dominance attribute, i.e. dominance 
rank, the simulation (2) that took distance of the asymptotic 
solution into consideration, and the simulation in which the 
selection process is controlled only by optimization criteria. 
Less effective were found the simulations (4) and (5) in which 
selection process was controlled only by optimization criteria, 
the simulation (6) that took into consideration the influence of 
dominance attribute, i.e. dominance count. Thus in the group 
of more effective strategies there were two simulations which 
use the combined fitness and one simulation which realizes the 
strategy of random combination of objective function without 


count 


dD Let’s notice that fulfillment of mentioned algorithm quality aspects can be considered as multi-objective task of algorithm optimization 


(optimizer). 
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Fig. 42. Specification of the sets of non-dominated solutions obtained during the performed genetic multi-objective optimization simulations 
of ship structure with respect to the structure weight f, and the surface area fy; black circles represent non-dominated solutions, 
red dots represent non-dominated solutions closest to the asymptotic one 
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taking into consideration domination attributes. On the basis of 
such superficial quality analysis can be proposed the statement 
that the combined fitness multi-objective optimization 
algorithm can be a more efficient strategy for multi-objective 
optimization of ship structure than scalarization of objective 
function strategy, without taking into consideration domination 
attributes in selection process. 

The conducted series of computer simulations confirmed 
efficiency of the developed computational algorithm and 
computer program for solving the formulated unified problem 
of topology and optimization of ship structure dimensions. As 
aresult of the calculations, were elaborated approximations 
of the set of Pareto-optimum solutions which in particular 
contain simulations of non-dominated solutions, from several 
to adozen or somewhat more in number. The obtained 
results do not allow to unequivocally determine superiority 
of any examined strategies of fitness function evaluation. For 
formulating more detailed quantitative conclusions further 
systematic statistical studies performed on much larger number 
of samples, are necessary. 

In the case of ship structure optimization the problem of 
efficiency assessment of elaborated algorithms is additionally 
complicated by the necessity of formulating many constraints 
and including them to the fitness function. In the task considered 
in this publication forty constraints were assumed. The 
constraint represented in optimization model in the form of 
proper components of penalty function, n, = 40, very strongly 
limit solution space available for searching; this way they 
also distort image of influence of the assumed fitness function 
computing strategies on algorithm convergence and the quality 
of attained solutions and therefore on the efficiency of the 
studied algorithms. 

From practical point of view it is interesting to notice which 
of the formulated constraints appears active in optimization 
process and which does not. In the first case to allocate large 
computational outlay to control them is justified. In the other 
case such outlay may appear useless from the point of view 
of effectiveness of optimization process and in some cases it 
is possible even to resign from controlling them. In the case 
of the optimization model used in this work the constraints in 
the form given by the inequality (21) concerning the required 
section moduli values of web frames in the three structural 
regions: side outboard region, bottom region and wet deck 
region, appeared active. In the case of the remaining constraints 
dealing first of all with the required thickness values of plates 
and dimensions of frames, the formulated constraints are 
satisfied with a large excess. 

The final conclusion can be formulated as follows: one 
cannot formulate a finite number of quantitative measures which 
allow putting in order the set of Pareto-front approximating sets 
in relation to the quality, and therefore one cannot formulate 
objective quality/efficiency measures of the proposed multi- 
objective optimization evolutionary algorithms. Thus one 
cannot prove objectively and unequivocally the supremacy of 
one of the algorithms proposed and discussed in this work or 
realization strategies of one of them. It depends on a potential 
user whether he / she would consider the presented concept 
interesting, elaborate its computer realization and finally verify 
its efficiency in his / her specific case. 


SUMMARY AND CONCLUSIONS 


The problem of minimization of weight and total surface 
area of the complete three-dimensional midship block-section 
of the high- speed catamaran hull was presented and discussed 
in detail. The strength criteria for checking ship structure were 


taken from the selected classification rules. The calculation 
tool for solving the formulated unified problem of the multi- 
objective optimization of topology and scantlings of the sea- 
going ship hull structure was developed with the accuracy 
typical for the preliminary design stage. 

The application of the genetic algorithm concept to solve the 
formulated optimization problem was presented. In the study it 
was proved that the genetic algorithm allows to include, in the 
multi-objective optimization model, a large number of design 
variables of real ship structure. The introducing of constraints 
related to strength, fabrication and standardization is not 
difficult and may cover a more representative set of criteria. 

The aggregation method was proved effective even in the 
case of the fixed values of the weight coefficients since in the 
case of the constrained problem the components of the penalty 
function introduce a random influence to the fitness function. 
The method is thus closer to that based on the random weight 
coefficients of the optimization criteria. 

This author has discussed crucial role of Pareto domination 
relation in process of evolution of feasible solutions for ship 
structure towards Pareto-front containing non-dominated 
variants of the ship structure. Using the concept of d omination 
in the set of feasible solutions this author has proposed his 
own definitions of the concept of domination rank as well 
as domination count, enabling this way to take into account 
relation between a feasible variant and other feasible variants. 
Basing on the ideas the author has proposed an evolutionary 
algorithm for solving the problem of topology-size multi- 
objective optimization of hull structure of sea-going ship, which 
uses, in selection process, the combined fitness function which 
allows taking into account, in selection process, the following 
items: (1) optimization criteria, (2) dominance attributes, 
(3) distance to the asymptotic solution as well as (4) penalty 
functions for violating assumed constraints. The computational 
program which makes it possible to perform ship structure 
multi-objective optimization with an accuracy appropriate for 
preliminary design stage, was elaborated. By using the tool 
and elaborated computational model of hull structure, series 
of computer simulations were conducted for the fast catamaran 
passenger-vehicle ferry of Auto Express 82 design. The results 
of the performed computations and subsequent discussion gave 
reasons for the statement that the elaborated algorithm may be 
considered an efficient tool for multi-objective optimization of 
ship structures in the preliminary design stage. 

It should be remembered that the developed multi-objective 
optimization algorithm is based on the random processes 
therefore the obtained results should be interpreted in the 
statistical sense. It means that the simulations and their results 
may appear sometimes not representative. a slight change of 
the developed models or control parameters may result in 
a different course of simulation and lead to different results. In 
this context further systematic studies on algorithm efficiency 
controlled by a particular component of combined fitness 
function, are necessary. 

Further systematic investigations of effectiveness of the 
proposed strategies, including repeated computations different 
to each other only by evolution history, aimed at statistical 
confirmation of the effectiveness, are deemed necessary. 
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ABSTRACT 


The article presents the results of the research project concerning tip vortex cavitation. This form of cavitation 
is very important in operation of many types of rotary hydraulic machines, including pumps, turbines and 
marine propellers. Tip vortex cavitation generates noise, vibration and erosion. It should be eliminated 
or significantly limited during the design of these types of machines. The objective of the project was to 
develop an accurate and reliable method for numerical prediction of tip vortex cavitation, which could 
serve this purpose. The project consisted of the laboratory experiments and numerical calculations. In 
the laboratory experiments tip vortex cavitation was generated behind a hydrofoil in the cavitation tunnel 
and the velocity field around the cavitating kernel was measured using the Particle Image Velocimetry 
method. Measurements were conducted in three cross-sections of the cavitating tip vortex for a number of 
angles of attack of the hydrofoil and for several values of the cavitation index. In the course of numerical 
calculations two commercial CFD codes were used: Fluent and CFX. Several available approaches to 
numerical modeling of tip vortex cavitation were applied and tested, attempting to reproduce the experimental 
conditions. The results of calculations were compared with the collected experimental data. The most 
promising computational approach was identified. 


Keywords: rotary hydraulic machinery; cavitation; numerical methods; experimental techniques 


INTRODUCTION 


One of the most important problems in design of different 
types of rotary hydraulic machines, such as marine propellers, 
water turbines and pumps is the as accurate as possible 
prediction of cavitation properties of these machines. Such 
a prediction enables elimination or significant limitation of 
the cavitation phenomena at the design stage. Among different 
types of cavitation, the tip vortex cavitation plays a very 
important role, because it is responsible for generation of 
intensive noise and vibration and also often causes erosion of 
the machinery elements. Development of accurate and reliable 
method for numerical prediction of tip vortex cavitation is the 
objective of research described below. 

The mechanism of formation of the cavitating tip vortex 
is shown schematically in Fig. 1. Combination of the inflow 
velocity to the hydrofoil and intensive rotation of the liquid 
around the vortex leads to the region of strongly reduced 
pressure, with minimum at the centre of the vortex just behind 
the tip of the foil. The complicated contradictory phenomena 
of concentration and dissipation of vorticity plays an important 
role in this process. Cavitation nuclei, i.e. micro-bubbles 
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naturally present in the liquid, are pushed by the pressure 
distribution into the centre of the vortex, where in sufficiently 
low pressure they undergo rapid growth, leading to formation 
of the cavitating kernel of the vortex. 


Fig. 1. Scheme of formation of the cavitating tip vortex 


The accurate determination of the pressure distribution in 
close vicinity of the tip vortex is the obvious pre-requisite for 
effective prediction of the tip vortex cavitation. The methods 
for calculation of the velocity and pressure around a tip vortex 
were the subject of earlier research published in [1, 2, 3, 6]. 

The continuation research, presented in this article, 
consists of the experimental and numerical parts. The purpose 
of the experimental part is to supply information about the 
velocity field around the cavitating tip vortex, necessary for 
development of the numerical methods and data concerning the 
geometry of the cavitating vortex for experimental verification 
of the numerical results. In the numerical part of the research 
project two commercial computer CFD codes were used and 
several available models of turbulence were tested. The detailed 
presentation of the experimental and numerical research is 
included in the following sections of the article. 


EXPERIMENTAL MEASUREMENTS OF 
THE VELOCITY FIELD AROUND THE 
CAVITATING TIP VORTEX 


The experiments were conducted in the cavitation tunnel 
of the Department of Energy and Industrial Apparatus of the 
Gdansk University of Technology, shown in Fig. 2. This tunnel 
has a rectangular measuring section having the dimensions 
3.0*0.35*0.45 meters. The maximum flow velocity in the tunnel 
is 6 meters per second. The velocity measurements around 
the cavitating tip vortex were performed using the Particle 
Image Velocimetry method and they were conducted by the 
team from the Warsaw University of Technology, Faculty of 
Civil Engineering, Mechanics and Petrochemistry. The PIV 
equipment set-up is also shown in Fig. 2. 


The hydrofoil model used in the experiments was specially 
designed on the basis of the typical contemporary marine 
propeller blade geometry, which was developed onto a plane 
surface (cf. Fig. 3). The span of the hydrofoil was selected 
equal to 225 mm in order to ensure that the cavitating tip vortex 
was located approximately in the centre line of the measuring 
section. The hydrofoil was manufactured of bronze and 
mounted vertically in the tunnel by means of the mechanism 
enabling accurate control of the angle of attack. 

The PIV measurements were conducted for a number of 
conditions, resulting from combinations of angle of attack and 
flow velocity. These conditions are listed in Table 1, together 
with the locations of three measurement planes. These planes 
were perpendicular to the tunnel axis and they were located 
50, 200 and 300 mm behind the tip of the hydrofoil, in order 
to visualize the process of development of the cavitating tip 
vortex. The combination of all conditions listed in Table 1 
produced 27 measurements of flow velocity (3 velocities * 
3 angles of attack * 3 measuring planes). During the experiments 
the static pressure in the tunnel was kept constant at 15 [kPa] 
and the variation in the cavitation index were obtained through 
changes in the flow velocity. The cavitation index o is defined 
according to the formula: 


-Pep 
l x72 (1) 
z 

where 

p — the static pressure in the tunnel 

p, — the critical vapour pressure 

p — the liquid density 

V — the flow velocity 


Tab. 1. Conditions for PIV measurements 


Distance 
behind the 
hydrofoil 


Cavitation 
indices 


Flow 
velocities 


tip [m/s] [-] 


4.32, 5.09, 5.87 
4.32, 5.09, 5.87 
4.32, 5.09, 5.87 


1.393, 1.003, 0.755 
1.393, 1.003, 0.755 
1.393, 1.003, 0.755 


4,8,12 
4,8, 12 


PIV method enables determination of the velocity vectors 
on the basis of measurement of displacement of particles 
between two correlated pictures created by the laser light 
(cf. Fig. 5) and registered consecutively in a short time interval 


Perspective 


Fig. 3. Design of the model hydrofoil 
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by the camera. In order to provide sufficient number of micro- 
particles in the flow, necessary for effective PIV measurements, 
the water in the tunnel was seeded with a silver metallic paint. 
The deflection of the light beam when crossing the window- 
water boundary was taken into account. The results of PIV 
measurements were stored in computer files and then filtered 
and re-calculated using specialized software. 

The scheme of the PIV measuring system is shown in Fig. 4. 
The distances of the camera K from the measuring section were 
kept constant b = 570 mm and c = 668 mm. The distance of 
the laser L from the measuring section was constant and equal 
to a= 440 mm. The distances of the measuring plane from the 
hydrofoil tip were varied: d = 50, 200 and 300. The system 
was arranged in such a way that the angle between the laser 
light plane and the camera axis was approximately equal to 45 
degrees in all measurements. 


2 i 
J; 
Y F 


H 


Fig. 4. Scheme of the measuring system 


L] 


Fig. 5 shows the laser light in the measuring plane 50 mm 
behind the hydrofoil tip during the PIV measurements. Apart 
from the PIV measurements the photographic registration of the 
cavitating tip vortex was performed. The complete description 
of the measurements may be found in [5]. 


Fig. 5. The PIV measurement plane behind the tip of the hydrofoil 


NUMERICAL PREDICTION OF THE 
CAVITATING TIP VORTEX 


The purpose ofthe first stage of calculations in the research 
project, described in this article, was to test and compare 
different CFD programs and different turbulence models 
from the point of view of their ability to predict accurately the 
geometry of the cavitating kernel of the tip vortex. One flow 
condition was selected for this comparison, namely hydrofoil 
angle of attack equal to 8 degrees and the velocity of flow 
equal to 5.2 m/s. The geometry of the computation domain 
taken into account in the CFD calculation is shown in Fig. 6 It 
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may be seen that quite a large section of the cavitation tunnel, 
especially in front of the hydrofoil, was modeled numerically. 
The inlet confuser was taken into account in the simulations in 
order to obtain the similar velocity distribution upstream the 
hydrofoil as in the measurements. Calculations were performed 
using two commercial CFD programs: Ansys/Fluent v12 and 
Ansys/CFX. In order to demonstrate the practical accuracy of 
numerical prediction of the cavitating tip vortex, the authors 
of calculations were not familiar with the experimental results 
beforehand. They prepared and performed the calculations 
using their best experience. 

The unstructured computational grid for Fluent is 
constructed of about 4.7 million of hexahedral elements. The 
mesh was created in Hexpress/Numeca mesh generator. The 
density ofelement distribution is increased near the tunnel walls 
in order to keep y=1. In view ofthe anticipated presence of the 
tip vortex the grid had also significantly increased density in the 
region behind the hydrofoil tip — see Fig. 7. In computations 
the MUSCL (Monotone Upstream-Centered Schemes for 
Conservation Laws) scheme was used and 4 turbulence models 
were tested: 

- standard k — s model 
- k-—eRNG 
- k-@SST 
- Reynolds Stress Model (RSM) 
Boundary conditions were set as follows: 
— at the inlet plane: 
- mass flow rate 437 kg/s (it corresponds to velocity at 
the test section ~5.2 m/s) 

- total temperature 283 K 

- turbulence intensity 1% 

- turbulent viscosity ratio 10 
— atthe outlet plane 

- static pressure 15 kPa 


ZX 


Y 


Fig. 6. The domain of flow taken into account in CFD calculations 


Fig. 7. Discretization of the hydrofoil and the plane behind it used in Fluent 


A mixed two-phase model was used in the calculations of 
the cavitating flow. Cavitation was determined on the basis of 
Rayleigh-Plesset equation in Zwart-Gerber-Belamri formulation 
[7]. The gaseous phase was treated as a compressible medium 
according to the perfect gas model. 

Calculations performed using the Ansys/CFX program 
were performed on the basis of unstructured grid having about 
9 million elements, including about 8.2 million tetrahedral 
elements and about 0.8 million prismatic elements in the 
boundary layer. Part of this grid located in the vicinity of the 
hydrofoil is shown in Fig. 8. The gaseous phase was treated as 
a incompressible medium in case of Ansys/CFX simulations. 


Fig. 8. Discretization of the hydrofoil and plane behind it used in CFX 


In the calculations by CFX the same turbulence models 
were used as in Fluent case. 

The comprehensive description of CFD calculations may 
be found in [4]. 


COMPARISON OF THE EXPERIMENTAL 
AND NUMERICAL RESULTS 


Comparison of the registered and calculated 
geometry of the cavitating tip vortex 


The photograph of the cavitating tip vortex in the analyzed 
flow condition (angle of attack 8 degrees, velocity of flow 5.2 


a 


pa 
“y 


300mm = 200mm 50 mm 


m/s) is shown in Fig. 9, while the results of corresponding CFD 
calculations are presented in Figs. 10, 11 and 12. The detailed 
analysis of the figures leads to the following observations: 

- there are significant differences between Fluent and CFX 
predictions, even when using the same turbulence model; 
these differences may be attributed to the markedly 
different structure in the computational grid between both 
programs, 

- similarly important differences may be found in predictions 
by each of the programs (Fluent or CFX) while using 
different turbulence models, 

- CFX results are much less dependent of the turbulence 
model than the Fluent ones, 

- the differences in performance between Fluent and CFX 
for the same turbulence models result most likely from 
the different structures of the computational grids and 
numerical schemes implemented in the solvers 

- k-e RNG turbulence model seems to be the best suited for 
prediction of tip vortex cavitation and it performs equally 
well both in Fluent and CFX, 

- standard k-s turbulence model in Fluent case gives the 
shortest cavitating zone what is highly related to the 
dissipative nature of this model 

- CFX results for standard k-e and SST turbulence model 
give surprisingly the same results 

- both programs predicted some sheet cavitation at the root of 

the hydrofoil which was not observed in the experiment. 


Fig. 9. Photograph of the cavitating tip vortex in the selected flow 
configuration 


Comparison of the measured and calculated 
velocity field near the cavitating vortex 


The measurements and calculations of the velocity field in 
the close vicinity of the cavitating tip vortex were performed 
for the same selected flow condition. The results are presented 
in the form of velocity vectors in the X-Y plane (i.e. plane 


300mm 200mm 


50 mm 


Fig. 10. Calculated cavitation - k-e standard turbulence model — Fluent (left) and CFX (right) 
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Fig. 11. Calculated cavitation — k-e RNG turbulence model — Fluent (left) and CFX (right) 


Fig. 12. Calculated cavitation — k-w SST turbulence model- Fluent (left) and CFX (right) 


perpendicular to the axis of the cavitating vortex) and they are 
grouped together according to the distance of this plane from 
the hydrofoil tip. Figs. 13, 14, 15, 16 show the results for 50 
mm distance, Figs 17, 18, 19, 20 — for 200 mm distance and 
Figs. 21, 22, 23 and 24 — for 300 mm distance. The results of 
measurements are presented (Figs. 13, 17 and 21) according to 
two different procedures of averaging (left and right section of 
the Figures) of the velocity field. The processing of experimental 
data in each flow conditions was done by PIV measurements of 
approximately 20 to 25 instantaneous velocity fields and their 
time averaging. Two final types of velocity fields were obtained 
for each measuring conditions by subtracting either space 
averaged velocity vector or only the horizontal component of 
this vector from the time averaged velocity field. The origin 
of the coordinating system for each PIV intersection distance 
from hydrofoil was in taken in the center of the cavitating tip 
vortex. The homogeneous high velocity area in the vicinity of 
the origin should be analyzed with respect of remarks described 
in the conclusions section of this article. 

Close inspection of the results for the plane located 50 
mm behind the hydrofoil tip leads to several interesting 
observations: 

- measurements indicate presence of two almost equally 
strong vortex kernels; such a situation may happen in 
specific flow cases, but in this case it is not confirmed by 
photographs and by calculations, 

- predictions by CFX look almost the same for all turbulence 
models, what corresponds to the results of calculation of 
cavitation described above, 

- predictions by Fluent depend strongly on the turbulence 
model, 

- predictions by Fluent and CFX agree quite well with each 
other for k-e RNG turbulence model, 

- Fluent results show existence of the wake downstream of 
the hydrofoil, such effect is less visible in prediction by 
CFX, what arises from the mesh structure. 
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Analysis of Figs 17 — 20, describing the situation at 200 mm 
behind the hydrofoil, may be summarized in the following way: 
- the measurements indicate the advancing process of 

concentration and merging of the two initially detected 

vortex kernels into one, 

- calculation by Fluent and CFX indicate the process of 
turbulent dissipation of vorticity and resulting weakening 
of the vortex, 

- this calculated process of dissipation of turbulent vorticity 
depends strongly on the model of turbulence; it seems to be 
overestimated in case of standard k-s and k-œ SST models, 
what leads to prediction of too short cavitating tip vortex 
(cf. Figs. 9 — 12), because the calculated intensity of the 
vortex falls too rapidly with distance from the hydrofoil, 

- the k-e RNG turbulence model produces very similar results 
both in Fluent and CFX, moreover, they agree reasonably 
with the results of measurements as far as the maximum 
values of velocity are concerned. 

- Fluent results still indicate existence of the wake downstream 
of the hydrofoil, its interaction with the vortex results in an 
asymmetric structure of the vortex, it is better shown for k-e 
RNG and SST models (less dissipative than standard k-e). 


The situation at 300 mm behind the hydrofoil tip, shown in 

Figs. 21 — 24, may be summarized in the following way: 

- measurements show a single strong vortex of a rather 
unnatural rectangular cross-section, 

- as far as the calculations are concerned, only the results 
obtained from Fluent with k-e RNG turbulence model 
produce maximum values of velocity comparable with that 
determined in the measurements, 

- all other calculations predict too intensive dissipation of 
vorticity, what is confirmed by the cavitation prediction 
shown in Figs. 10 — 12 

- the interaction of the tip vortex with the wake is still present 
in Fluent results. 
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Fig. 15. Calculated velocity field in the plane X-Y at the distance 50 mm behind the hydrofoil tip — k-e RNG turbulence model — Fluent (left), CFX (right) 
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Fig. 16. Calculated velocity field in the plane X-Y at the distance 50 mm behind the hydrofoil tip — k-w SST turbulence model — Fluent (left), CFX (right) 
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Fig. 17. Measured velocity field in plane X-Y at the distance 200 mm behind the hydrofoil tip 


ooo S = Sane 
MVAMOHWDONM AN 


008 O1 O12 014 O16 018 12. 0. 
X [m] X [m] 


Fig. 18. Calculated velocity field in the plane X-Y at the distance 200 mm behind the hydrofoil tip — k-e standard turbulence model — Fluent (left), CFX (right) 
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Fig. 20. Calculated velocity field in the plane X-Y at the distance 200 mm behind the hydrofoil tip — k-w SST turbulence model — Fluent (left), CFX (right) 
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CONCLUSIONS 


The detailed analysis of the results of measurements, 
observations and CFD calculations presented in this article 
leads to the following conclusions: 

- the contemporary commercial CFD codes are generally 
capable of predicting the geometry of the cavitating tip 
vortices generated by hydrofoils reasonably well, 

- the accuracy of CFD prediction of the geometry of the 
cavitating tip vortex depends strongly on the turbulence 
models and on the grid structure, hence only the grids 
constructed especially for vortex-dominated flows should 
be used, together with turbulence models especially suited 
for modeling of such flows [1, 2, 3, 6], 

- inview of the above presented results the turbulence model 
k-e RNG seems to be best suited for CFD prediction of tip 
vortex cavitation, despite the fact that according to [3] the 
k-@ SST performed best in prediction of non-cavitating 
vortex flows, 

- the measurements of the velocity field in the vicinity of 
the cavitating tip vortex by means of PIV methods seem 
to be a difficult and challenging task, especially due to the 
following reasons: 

o unsteady oscillations of the cavitating kernel of the 

vortex, 

o uncontrolled content of cavitation nuclei (i.e. gas and 
vapour filled micro-bubbles) carried by the flowing 
water, 

o shading of part of the measuring plane by the cavitating 
kernel, 

o difficulties in homogeneous PIV seeding due to the large 
volume of fluid inside the cavitation tunnel and complex 
flow conditions near cavitating tip vortex (centrifugal 
force acting on seeding particles) 

- due to the above listed reasons the velocity distribution 
predicted numerically at the following sections is 
significantly different from the measurements. 
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ABSTRACT 


Free vibration of grid stiffened circular cylindrical shells is investigated based on the first Love's 

approximation theory using Galerkin method. Full free edges are considered for boundary conditions. An 

equivalent stiffness model (ESM) is used to develop the analytical solution of the grid stiffened circular 

cylindrical shell. The effect of helical stiffeners orientation and some of the geometric parameters of the 

structure have been shown. The accuracy of the analysis has been examined by comparing results with 
those available in the literature and finite element approach. 


Keywords: Free vibration; Grid Structure; Cylindrical Shells; Stiffener; Galerkin Method 


INTRODUCTION 


Grid stiffened cylindrical shells are applicable in many 
industries such as aerospace industries. Stiffened cylindrical 
shells play a big role in aerospace industries in fuselage and 
fuel tank applications. This has resulted in an extensive research 
work in the field of cylinders with stiffening structures [1-8]. 

Grid stiffened cylindrical shells are widely used in 
engineering fields. These structures are subjected to external 
dynamic loads. These external dynamic loads can cause the 
undesirable resonance and it can lead to fatigue. Moreover, 
dynamic characteristics must be used on design of structure 
because only vibration (not fatigue) could severely damage 
the sensitive equipment in airplanes, launch vehicles and etc. 
Therefore, it is essential to understand the dynamic behavior 
of these structures. Theoretical methods of analyzing the 
grid stiffened structures are classified into two main types, 
depending upon whether the stiffeners are treated by averaging 
their properties over the shell surface to conventional materials 
or by considering them as discrete elements. The first method, 
so-called smeared stiffener theory, is particularly applicable 
only when large numbers of stiffeners are closely and evenly 
spaced. The second method, so-called discrete stiffener theory, 
is more general as it can accommodate any stiffener distribution. 
Numerous researches have been developed to study the 
vibrational behavior of stiffened cylindrical shells. Mustafa and 
Ali [9] predicted natural frequencies for the stiffened cylindrical 
shells using the Rayleigh-Ritz procedure. In this procedure they 
used only one term in assuming the displacement functions 
satisfying the simply supported boundary condition. One-term 


approximation is sufficient for the analysis of the cylindrical 
shells with simply supported boundary condition. However, it 
can lead to much error to obtain the exact solution of stiffened 
shells with any other boundary conditions. Yang and Zhou [10] 
presented the transfer function method to analyze the ring- 
stiffened shell. Lee and Kim [11, 12] investigated the effect 
of rotation speeds and boundary conditions on the frequencies 
for the orthogonally stiffened composite cylindrical shells 
treating the materials of stiffeners as equivalent isotropic. The 
mentioned papers were, however, limited to the shells with the 
uniform dimensional and evenly spaced stiffeners. In fact, non- 
uniform dimensional and unevenly spaced stiffeners are used 
much more in structural reinforcements. Wang et al. [13] solved 
the free vibration problem for the isotropic cylindrical shells 
with varying ring-stiffener distribution using the extended Ritz 
method. Egle and Sewall [14], in different boundary conditions, 
have analyzed the effect of stiffeners on natural frequencies 
of stiffened cylindrical shells. In this research, stiffeners are 
considered as discrete elements, energy method and Hamilton 
principle are used to obtain equations of motion. Rinehart and 
Wang [15] have studied the changes in natural frequencies 
of cylindrical shells affected by stringers stiffeners, based 
on Vlasov thin walled beam theory. By both, considering the 
stiffeners as discrete elements and using energy method, they 
have obtained the equations of motion. 

Unlike the previous study, the vibrations of stiffened 
cylindrical shells with grid structure under full free boundary 
conditions are analyzed in this paper. The stiffness matrix of 
the whole structure is determined by stiffness matrix of grid 
structures. Then, equilibrium equations are considered based 
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on the classical shell theory. Strain-displacement relations are 
written based on the first approximation of Love theory and 
then by replacement in stress-strain relations, equilibrium 
equations based on displacement parameters are obtained. After 
simplification of equilibrium equations, the shell frequency 
equation is obtained by using Galerkin method. Finally, 
according to the following assumptions, the effect of stiffener 
geometry and mass of grid stiffened cylindrical shell on natural 
frequencies is presented: 

1. The thickness of ribs is small compared with length of the 
ribs. So the transverse strain of stiffeners is much smaller 
than that of longitudinal strain and can be negligible. 

2. The strain is uniform across the cross-sectional area of the 
stiffeners. Hence, a uniform stress distribution is assumed 
across the cross-sectional area of the stiffeners. 

3. The load on the stiffener/shell is transferred through shear 
forces between the stiffeners and shell. 
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Fig. 1. Grid stiffener and shell design parameters 


The above matrix elements are the functions of the strain 
and curvature parameters of the shell. We denote these 
stiffness parameters by As, Bs and Ds which corresponds to the 
extensional, coupling and bending matrices, respectively. Here, 
the superscript s stands for stiffener. By the matrix summation 
of force and moment, coming from the stiffener Grid structure 
and cylindrical shell, the force and the moment of the whole 
structure will be obtained. 


Me N° +N" 
M| |M"+ M" 


Ns and Ms are the force and moment contribution of the 
shell, respectively. Forces and moments, which affect the shell, 


(3) 
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ANALYTICAL METHOD 


Equivalent Stiffness Matrix 
of the Stiffener and Shell 


To calculate the natural frequencies of the whole structures, 
first, the stiffness matrix of grid structures must be determined. 
This matrix is composed of three separated matrices As, Bs and 
Ds which can be explained like below: 


NP [AS B5 |fe? 
M| |B ps|| k 
Based on geometrical variables of grid structures which 


presented in Fig. 1, stiffness matrix of grid structures can be 
obtained as below [16]: 


(1) 


3 2 
A Sg 
a a 
2 3 
cst s +L 0 r _ 
b b x 
2 ge 
0 0 s“ct : 
a Vx 
e415 a k, (2) 
t cs — 0 : 
2a 2a ko 
2 2 
s’?’— ($+1)— 0 Kox 
2b 2b 
0 0 ge 
2a 


relate to the occurred strain by stiffness matrix A*", Bs and D+, 
Based on shell mechanical properties, these matrices are: 


m lv 0 
Av=— |v 1 0 | BS=[0] 
“ l-v 
2 (4) 
ae lI v 0 
—— v 0 | i j=1,2,6 
* dev) 
0 0 — 
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Substituting the force and moment expressions for the 
stiffener network from Eq. 2, and the force and moment 
expressions for the shell from Eq. 4, the total structure 
constitutive equation given by Eq. 5 results: 


N]_ A+A? B+B” |e 
M| [B'"+B* D*+D* || k a 
Formulation 


The cylindrical shell under consideration is with constant 
thickness t, radius R and length L. The reference surface of the 
shell is taken to be at its middle surface where an orthogonal 
coordinate system (x, 0, z) is fixed. As shown in Fig. 2, the 


x axis is taken in the axial direction of the shell, where the 8 and 
z axises are in the circumferential and radial directions of the 
shell, respectively. The displacements of the shell are defined 
by u, v, w in the x, 0, z directions respectively. 

The equations of motion for a cylindrical shell can be 
written by the Love theory as: 
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Fig. 2. Co-ordinate system and circumferential modal shape [17] 
Boundary Conditions 


Due to the satisfaction of the boundary conditions, the 
displacement u, v and w can be explained as double Fourier 
series [19]: 


u(x,0,t)= r> Aun Au cos(n0) cos(t) 
x 


v(x,0,t)= >>, Bw? (x)sin(n8)cos(@t) (7) 


W(x,9,t) = >> Can P(X) Cos(nO ) cos(wt) 


In recent equations, A,,,, B,,, and Can are coefficients 
of natural modes’ shape, which obtained from solving free 
vibration. For solving free vibrations, T (t) = e®mn! is considered 
as a function of time, m is the number of axial half-wavelength, 
n is the number of circumferential half-wavelength and @,,,, 
is the natural frequency in mode of mn. To satisfy boundary 
conditions, axial and circumferential functions are explained 


as below: 
AX A Xx 
z |+a,cos| — 
L L 


(x)= a, cosh 


h 0) =sin(n0 ) 


In above equations a; are constant coefficients which 
determined according to boundary conditions. i, is the root of 
non-linear equations and o,, is the dependant parameter on An 
which obtained according to boundary conditions. Free-Free 
supported conditions can be defined as below: 


Hox) _ F(x) _ 
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Free Vibration Analysis 


Using Galerkin method and substituting Eq. 7 into Eq. 6, 
it can be written as: 
Cy Cy G3 JA j 
Ban f= 0 a 
Cs, Cx €33 | [Dann 
where: 
C,(i, j= 1, 2,3) — the parameters from the L; after they are 
operated with the x and 0. 


For non-trivial solutions, one sets the determinant of the 
characteristic matrix in Eq. 9 to zero: 


det(C,) = 0 (i, j = 1, 2, 3) (10) 
So the frequency equation can be obtained as: 
Bas + Bot + Bm? + B, = 0 (11) 


where: 
B,G,j = 1, 2,3) — the coefficients of Eq. 10. 


Solving Eq. 11, one obtains three positive roots and three 
negative roots. The three positive roots are the angular natural 
frequencies of the cylindrical shell in the axial, circumferential 
and radial directions. The lowest of the three positive roots 
represents the flexural vibration, and the other two are in-plane 
vibrations. 


NUMERICAL RESULTS AND DISCUSSION 


Numerical implementation of the present analysis was 
performed using general-purpose computation package 
MATLAB. To check the accuracy of the present analysis, 
the results obtained are compared with those in the [17]. 
A comparison of the values of the frequency of a free vibrating 
cylindrical shell with the F-F boundary conditions is given in 
Table 1. 


Tab. 1. Comparison of values of the natural frequency œ, for a cylindrical 
shell with F-F supported boundary conditions 
m = 1; R/h = 374; L/R = 2.63; h = 0.6477 mm; 


E = 70 GPA; v = 0.3; p = 2700 Kg/m3 
Reference [17] 


Present | Difference (%) 


m X 


. [À 
-asin i 


) (=u, v, w) (8.a) 


In the parameter, E is Young’s modulus of elasticity, v is 
the Poisson ratio, p is the density, R is the radius and o is the 
frequency. The comparisons are carried out for the parameter of 
L/R = 2.63 and for the cases of R/h = 374 and t = 6.477x10+ m. 
Using the method outlined earlier, numerical results are 
obtained for the six model of grid-stiffened cylinder with four 
families of ribs and geometric and mechanical parameters given 
in Tables 2 and 3. 
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Tab. 2. Dimensions of the considered structures (mm, kg) 


m 


stiffener 


19.794] 3.450 


5 
5 
5 
> 
5 


20.944] 4.430 |16.344 300}640/45s| 6 
22.505) 6.161 |16.344| 5 |300/653|/60s) 6 


Tab. 3. Mechanical properties of the models (Aluminum) 


E (GPa) v p kg/m? 
70 0.3 2700 
The obtained results are compared to the solution of the 


same problem when the considered structures are modeled with 
a finite element analysis package ANSYS. One of the generated 


Fig. 3. The mesh generated for a grid stiffened cylindrical shell using 
ANSYS 


Based on a technical report [18], natural frequencies of 
a grid-stiffened shell can be approximated by equivalent of 
thickness. In this method, a conventional shell with the same 
weight as grid part is predicted. An equivalent shell will be 
produced by the replacement of this shell with grid part. The 
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Fig. 4. Values of natural frequency for Model (1) with Free-Free boundary 
conditions in three methods. A FEM; e Present Method; 0 ETM (Equivalent 
Thickness Method); m = 1 
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obtained frequencies of this equivalent shell can be a proper 
approximation for natural frequencies of the main structure. In 
Figs. 4-9, a comparison between the results obtained by ETM 
(Equivalent Thickness Method), FEM and present method is 
presented. These figures have shown that the present method 


gives more exact answers than ETM. 
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Fig. 5. Values of natural frequency for Model (2) with Free-Free boundary 
conditions in three methods. A FEM; è Present Method; o ETM (Equivalent 
Thickness Method); m = 1 
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Fig. 6. Values of natural frequency for Model (3) with Free-Free boundary 
conditions in three methods. A FEM; @ Present Method; o ETM (Equivalent 
Thickness Method); m = 1 
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Fig. 7. Values of natural frequency for Model (4) with Free-Free boundary 
conditions in three methods. A FEM; @ Present Method; o ETM (Equivalent 
Thickness Method); m = 1 


Fig. 10 shows the changes relative to natural frequencies 
of six predicted models, according to the angle of helical 
ribs. As shown in this figure, when the angle of helical rib 
increases, the amount of natural frequencies of structures will 
be increased. 
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Fig. 8. Values of natural frequency for Model (5) with Free-Free boundary 
conditions in three methods. A FEM; è Present Method; o ETM (Equivalent 
Thickness Method); m = 1 
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Fig. 9. Values of natural frequency for Model (6) with Free-Free boundary 
conditions in three methods. A FEM; e Present Method; 0 ETM (Equivalent 
Thickness Method); m = 1 
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Fig. 10. Comparison of natural frequencies of three predicted models 
according to the angle of helical ribs. © œ = 30 (Model 1); A ọ = 45 (Model 
5); o g = 60 (Model 6); m = 1 


Note: In this comparison Model 5 is smaller length than the 
remaining ones (about 2%). According to difference between 
configurations of models, the resulted value from FEM modeling 
is the minimum height discrepancy can be obtained. 


CONCLUSION 


The article has presented the analysis of grid stiffened 
cylindrical shell using Galerkin method. Comparison of the 
results by the present method and numerical finite element 
method was carried out. The six finite element models for 
grid stiffened cylindrical shell were created. The shell was 
fully free-free at both ends. The second six natural frequencies 
were obtained with the ANSYS. These results were compared 


with the present method and the agreement between them was 
good. Based on comparisons of the mentioned method, it is 
concluded that the present method is more convenient, more 
effective and more accurate. 
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torque of mechanical losses in the pump used 
in a hydrostatic drive 
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ABSTRACT 


The paper presents theoretical and mathematical models of the torque of mechanical 
losses in the pump with theoretical (constant) capacity qp, per one shaft revolution (with 
constant theoretical working volume Vp) and geometrical (variable) capacity qp,., per one 
shaft revolution (with variable volume V>,,). The models may be used in the laboratory 
and simulation investigations of the pump energy efficiency and the hydrostatic drive 


efficiency. 


Key words: hydrostatic drive; displacement pump; energy efficiency 


INTRODUCTION 


The paper is a continuation of the work presented in 
references [1 + 18], aimed at creating a method of evaluation of 
the losses and energy efficiency of hydrostatic drives as well as 
the used displacement machines (pump and hydraulic motors). 
The method is based on mathematical models of energy losses 
in the pumps, in hydraulic motors and in other elements of 
a hydrostatic drive system. 

The description of pump losses and energy efficiency is 
based on the diagram of power increase in the drive system 
opposite to the direction of power flow, replacing the Sankey 
diagram of power decrease in the direction of power flow 
[18]. The Sankey diagram of decrease (division) of power 
in a drive system in the direction of power flow is the main 
reason of incorrect evaluation of the energy losses, a. o. in the 
displacement pumps and hydraulic motors used in hydrostatic 
drive systems. 

During the operation of a hydrostatic drive system, the 
energy losses force the increase of power in the system 
— from useful power Pm. required by the hydraulic motor 
driven machine to the power Pp, consumed by the pump on 
the pump shaft. 

In the description of power flow in a drive system and 
the powers of energy losses connected with the flow, the 
notions: „power decrease”, „power division”, „power loss” 
should not be used. 

The notion associated with the power of energy losses 
in a drive system is ,,increase of power”. 

Figure 1 presents a diagram of power increase in 
a displacement pump opposite to the direction of power flow, 
which replaces the Sankey diagram of power decrease in the 
direction of power flow. 
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The aim of the paper is to present the theoretical and 
mathematical models of mechanical losses in the pump 
„Working chambers — shaft” assembly. Pump is a displacement 
machine with theoretical (constant) capacity qp, per one shaft 
revolution (with constant theoretical working volume V>,) 
or with geometrical (variable) capacity qp,, per one shaft 
revolution (with variable geometrical working volume Vp). 

The models may be used in the laboratory and simulation 
investigations of the pump mechanical losses, allowing to 
evaluate the pump energy efficiency and the hydrostatic drive 
efficiency . 


THEORETICAL MODELS OF THE TORQUE 
Mpm OF MECHANICAL LOSSES IN THE 
PUMP „WORKING CHAMBERS - SHAFT” 
ASSEMBLY 


The pump shaft torque M, (required by the pump of its 
driving motor) must be greater than the torque M,, indicated 
in the pump working chambers because of the necessity of 
balancing also the torque M,,, of mechanical losses in the 
working chambers — shaft” assembly. The assembly forms 
the working chambers and changes their capacity and also 
connects the working chambers with the shaft. Therefore, 
the torque M, required on the pump shaft is a sum of the torque 
M; indicated in the working chambers and the torque Mpm of 
mechanical losses in the pump ,,working chambers — shaft” 
assembly: 


Mp = Mpi gs Mom (1) 
Torque M,,, of mechanical losses in a pump with variable 


capacity qp,, per one shaft revolution is, at the maximum value 
Of pey 1-€. Fpey = Gp: (With bp = gp,,/qp, =1), equal to the torque 


viscosity » of the working fluid 
inlet pressure pp, =0 
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Fig. 1. Diagram of power increase in a displacement pump opposite to the direction of power flow, 
replacing the Sankey diagram of power decrease in the direction of power flow 


of mechanical losses in that pump working as a pump with 
constant capacity qp per one shaft revolution. The theoretical 
and mathematical models describing the torque Mpm of 
mechanical losses ina pump with variable capacity qp,, per one 
shaft revolution may be based on models of Mp, describing the 
torque of mechanical losses in a pump with constant capacity qp, 
per one shaft revolution (with b, = 1). Considering the models 
describing the torque of pump mechanical losses, we assume, 
that the pump is driven with practically constant rotation speed 
np and the decrease of shaft speed (decrease of the pump driving 
motor speed as a result of the increase of torque M, loading 
the motor shaft) to a value np < npo (Np) — rotational speed of 
unload pump driving motor) is negligible from the point of 
view of the impact of speed np on the value of torque Mpm of 
mechanical losses. 

Torque M,,, of mechanical losses in the pump is mainly 
an effect of friction forces between elements of the „working 
chambers — shaft” assembly, depending, among other, on 
the torque M, indicated in the working chambers — M}; = 
= Gpgy APpi/2 11 = bp qr: App/2 11. 

Friction forces between elements of the „working 
chambers — shaft” assembly are, to some extent, also an 
effect of the load on those elements of the inertia forces 
from rotational and reciprocating motion and depend 
on the pump capacity qr, per one shaft revolution 
(bp coefficient). 

In the piston (axial or radial) pumps with casing 
(crankcase) filled with liquid, friction forces also occur 
between elements of the „working chambers — shaft” 
assembly and the liquid and depend on the liquid viscosity v. 

The value of torque Memjap).tpv, Of mechanical losses in 
the pump ,,working chambers — shaft” assembly, loaded with 
indicated increase App; of pressure in the working chambers, 
in the pump operating at the capacity qp,, = b, qp,per one shaft 
revolution and discharging the working liquid with (constant) 
reference viscosity v,, can be described as a sum of torque 
Mbmlappicbp.vn Of mechanical losses in the unloaded pump (torque 
of the losses when the indicated increase App; of pressure in 
the pump working chambers is equal to zero — App; = 0) and 
increase Mpmippi.tp.vn Of torque of mechanical losses, the increase 
being an effect of loading the assembly structure elements with 


torque M,, indicated in the pump working chambers (torque 
Mp; appearing when the indicated increase App; of pressure in 
the pump working chambers is greater than zero — App;> 0): 


(2) 


Torque M,,; indicated in the pump working chambers is 
proportional to the increase App, of pressure in the chambers 
and to the active volume of the chambers created during 
one pump shaft revolution, which is equal to the theoretical 
capacity qp, per one shaft revolution in a pump with constant 
capacity per one shaft revolution or to the geometrical capacity 
pe = bp gp, per one shaft revolution in a pump with variable 
capacity per one shaft revolution. 

Therefore, the „working chambers — shaft” assembly 
structure elements are loaded: 

— inapump with theoretical (constant) capacity qp, per one 
shaft revolution — with indicated torque: 


_ Wp APpi 
My = 
211 
— ina pump with geometrical (variable) capacity qp,, per one 
shaft revolution — with indicated torque: 


_ IrevAPpi _ bpqpiAPpi 


Pi 0m7 20 


which, in effect, causes a differentiated intensity ofthe increase 
Mbmlappicbp.vn Of the torque of mechanical losses, determined, with 
different values of coefficient bp = qp,,/qp,, as a function of the 
increase App; of pressure in the pump working chambers. 

In the theoretical and mathematical models describing 
the torque Mpmjappibpvn Of mechanical losses a hypothesis 
is assumed, that the increase Mpmippibpvn Of the torque of 
mechanical losses in the pump is proportional to the torque 
Mx indicated in its working chambers (Fig. 2 and 5). 

The impact of inertia forces of the „working chambers 
— shaft” assembly elements, performing the rotational and 
reciprocating motion in the pump, on the torque Mpm of 
mechanical losses can be presented, under the assumption 
that rotational speed n, of the pump driving motor 
changes only in a small range, as a function of capacity qp,, 


Mom|app, bin Mo mlApp;=0,bp.v, T AM m|ap,; Drv 


POLISH MARITIME RESEARCH, No 4/2011 


29 


M | APp; i Npo=cte, ps (bp =1 ) Va 
Pm|b, =1] 
A Ap App =0 Ap 
AD; APs = APpi 
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APpi=P, 
AMp bp =1 
Ya 
| \ App; =P, 
| App; =0 Mom bp =1 
Mem bp =1 Vr 
Yn 
| ! — 
0 Pn AP»; 
Fig. 2. Torque Mem\spp.b)=1.vn Of mechanical losses in the pump with constant capacity qp, per one shaft revolution (bp =1), 
with working liquid reference viscosity v„ as a function of the indicated increase App; of pressure 
in the pump working chambers — graphical interpretation of theoretical model (2) 
|Ap.,=0 | Mpo=Cte, App: =0, qa (bp =1) 
Mom | bo =1 
? 
| App, =0 
Mem |bp =1 
[Pn 


0 Vmin 


Vr 


Fig. 3. Torque Mrmsppi=0.bp-1» of mechanical losses in a piston (axial or radial) pump with crankcase filled with liquid and with constant capacity qp, per one 
shaft revolution (bp = 1), at the indicated increase App; = 0 of pressure in the pump working chambers, as a function of the ratio of viscosity v to reference 
viscosity v, — v/v,— graphical interpretation of theoretical model (3); torque Mpm\4pp;-0,b=1, Of mechanical losses in the pump without the crankcase filled with 
liquid is practically independent of the liquid viscosity v and is determined at the liquid reference viscosity v, 


(bp coefficient) per one shaft revolution of a variable 
capacity pump. Inertia forces do not depend on the value 
of increase Ap,; of pressure in the working chambers, 
therefore their impact on the torque M,,, of mechanical 
losses in the pump may be included in the evaluation of the 
torque Memippi-0,bp.vn Of mechanical losses determined at the 
increase App; = 0 (Fig. 5). 

The impact of the friction forces between the „working 
chambers — shaft” assembly elements and the liquid in the 
casing (crankcase) of the piston pump on the torque Mpm 
of mechanical losses in the pump can be presented, under 
the assumption that speed n, changes in a small range, as 
a relation of Mp,, to the viscosity v and to the capacity qp,, 
(b, coefficient) per one shaft revolution (Fig. 3, 4, 6, 7). 
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It is assumed, that the impact of liquid viscosity v 
on the friction forces between the „working chambers 
— shaft” elements and the liquid in the piston pump casing 
(crankcase), and in effect on the torque Mpn of mechanical 
losses in the pump, can be evaluated at one level of the 
increase App; of pressure indicated in the working chambers, 
e.g. at the increase App; = 0 (Fig. 3, 6). This assumption is 
connected with a simplification assuming that there is no 
significant impact of the increase App; of pressure on the 
liquid viscosity v and with assuming in the model describing 
the torque Mpm of mechanical losses the liquid viscosity v 
determined in the pump inlet conduit [at pressure pp, equal 
to zero (at liquid absolute pressure equal to atmospheric 
pressure) ]. 


The impact of inertia forces of structure elements performing 
the rotational or reciprocating motion in the pump and also the 
impact of liquid viscosity v on the torque Mpm of mechanical 
losses in the pump is then described in the model of the torque 
Mpmjappi-0.bpv Of those losses in an unloaded pump (at App; = 0) 
supplied with working liquid of changing viscosity v. 

The proposed theoretical models describing the torque 
Mbmiappi-0.bp.v of mechanical losses in an unloaded pump (at 
the indicated increase Ap,; = 0 of pressure in the working 
chambers) and at changing working liquid viscosity v [the 
impact of liquid viscosity v occurs in the piston pumps with 
liquid filling the casing (crankcase)] have the form: 

e inthe pump with theoretical (constant) capacity qp, (bp = 1) 

per one shaft revolution (Fig. 3): 


v 


Pm|App; =0, bp =1, Vp an, 
n 


M M 


(3) 


Pm|App; =0,bp =v 


e in the pump with geometrical (variable) capacity qp,, 
(Gp, = bp qp,) per one shaft revolution (Fig. 6): 


Momlapp, =0, bp, v = (4) 
aym 
=(M AM . 
= Pm|App; =0, bp =0, vat Pm|App; =0, bp, 2) ae, 
n 
where: 
Pm|App;=0, bp. Va 
= Momlapp;=0. bp, Yn Momlapp, =0, bp=0, v, = (5) 


= (Momlap,, =0, bp=l, vp MPmjappi =0, bp=0, “) bp 


Exponent a,,, in expressions (3) and (4) describes the 
impact of the ratio v/v, of working liquid v to reference 
viscosity v, = 35mm?s- on the value of torque of mechanical 
losses in a piston displacement machine with liquid filling 


App; i Npp=Cte, Gpy (bp =1 ) Vmin, Yn, Ymox 
Mom) be =! 
v 


App,=0 
Mom bp =1 
Vmax 


APp;=0 
bp =1 
Yn 


Mom 


App,=0 \ 
Mpm bp =1 


Vmin 


the casing (crankcase) (in the pump and in the hydraulic 
motor). 

The increase Mpmap,i.bpv Of the torque of mechanical 
losses in the pump, due to the load of the assembly elements 
with the indicated torque M), resulting from the indicated 
increase App, of pressure in the pump working chambers, 
is independent of the inertia forces of elements performing 
the rotational or reciprocating motion in the pump. It is also 
practically independent of the working liquid viscosity v; 
therefore, it may be determined at one viscosity value, e.g. 
at the liquid reference viscosity v, (Fig. 4, 7). 

The proposed theoretical models describing the increase 
Mbmippi,bpv Of the torque of mechanical losses in the pump, 
resulting from the indicated increase Ap,, of pressure in the 
working chambers, have the form: 

e inthe pump with theoretical (constant) capacity qp, (bp = 1) 

per one shaft revolution (Fig. 4): 


AM = AM 


Pm|App;.bp=l.v Pm|App;.bp=Lvn 


E Mjane bp1,¥, > Momlapp; =0,bp=l,v, (6) 


App; 


n 
e in the pump with geometrical (variable) capacity qp,, 
(drev = bp qr) per one shaft revolution (Fig. 7): 


AM AM 


= (Momlapp, =p,-bp=Lv, _ Momlapp: =0,bp=Lv, ) 


Pm|App;.bp.v Pm|App;.bp.V, 


E Momlapp, „Dp: Vn 7 Momlapp, =0, bp Vn 


=(M (7) 


+ 
Pm|App; =Pn- bp =l, v, 


n 


App; 


n 


-M 


Pm|App;=0. bp =1, v dbp 


Ym 
APp; APp; :0 P App 
Mpmlbp =1 =Mpm bp =1 7, +AMpm bp =I 
Vmax Va Va 


/ APp=Ph 
/ Ap App =0 Ap Membe =I 
/ Pp; Pi~ Appi 
/ Mpm|bp =1 =Mpm|bp=1 +AMpm bp =1 Pn 
/ Va Vn Yn 
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0 
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Fig. 4. Torque Mpm\dpp;bp=1v of mechanical losses in a piston (axial or radial ) pump with crankcase filled with liquid and with constant capacity qp, per one 
shaft revolution (bp = 1), as a function of the indicated increase App, of pressure in the pump working chambers — graphical interpretation of theoretical 
models (2) and (8); liquid viscosity Vinin Vn ANA Ving Torque Mpmappi.bp=1.v of mechanical losses in the pump without the crankcase filled with liquid is 
practically independent of the liquid viscosity v and is determined at the liquid reference viscosity v, 
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In effect, the proposed theoretical models describing 
the torque M,,, of mechanical losses in the pump take the 
forms: 

e inthe pump with theoretical (constant) capacity qp, (bp = 1) 

per one shaft revolution (Fig. 4): 


Momlapp:.bp =l,v— Momlapp;=0,bp =Lv, E + 
Vn 
(8) 
aym 
i AMbmjappi „bp =l,vn Momlape: =0, bp =l, vy ae = 
n 
APp; 
Pi 
= (Momlape, =p, -bp =l, Va = Momlapp; =0,bp=l.v, p 
n 


e in the pump with geometrical (variable) capacity qp,, 


(drev = bp qp,) per one shaft revolution (Fig.7): 
Momlapp: „Dp V = (9) 
aym 
=(Momapp,=0.bp =0,v, + AMpmļAppi=0,bp,v, ) + 
n 

+ AMpmlap»; „bp, Vn 
where: 

AM 


Pm|App;=0.bp; Vn = 


= (Mbmjapp: =D, bpl vy Mo mlapp, =0, bp =0, a) bp (10) 


AM 


Pm|App;-bp. Vn 
AP pi 


n 


(M 


Pm|App; =Pn> bp =l, Va Mpmjappi =0, bp =l, v Jbp a 1) 


MATHEMATICAL MODELS OF THE 
TORQUE OF MECHANICAL LOSSES 


In the mathematical models describing the torque Mpm 
of mechanical losses in the pump, coefficients k; of losses 
are used relating (comparing) the components describing 
the torque M,,, of losses in theoretical models to the pump 
theoretical torque M,,. The pump theoretical torque Mp, is 
also a reference value used in the description of the torque Mp, 
indicated in the pump working chambers: 

— theoretical torque: 


= AptPn 
211 


of the pump, with theoretical (constant) capacity qp, per one 
shaft revolution (bp = 1), is determined with the increase 
App of pressure in the pump equal to the system nominal 
pressure p, — App = p,, and with the assumption that there 
are no pressure and mechanical losses in the pump, 

— indicated torque: 


Mp, = qptAPpi — FptPn APpi =Mp, App; 


~ 20 2p, Pn 
in working chambers of the pump with theoretical (constant) 
capacity qp, per one shaft revolution (b, = 1) is determined 
with the indicated increase App; of pressure in the working 
chambers, 
— indicated torque: 


Pt 


M- qdpPevAPri _ bpgpAPpi _ 
Pi =~ -zm 
21 211 
= ptPn b Appi =M.„b APpi 
a p T 


in working chambers of the pump with geometrical 
(variable) capacity qp,,= bp qp; per one shaft revolution is 
determined with the indicated increase App; of pressure in 
the working chambers. 


Mpmlt eead Mog=Cle, Qag =0 (b=0), qeg (bp), drg =q (b,=1), Va 
mlb, 
Me 
AD»; App, =0 A 
Mom Dp 24 =Mpm oy +HAMpm er 
Vr Pa z 


Moni Z0 Meny i ° bo AMen EP 


Ya Va 
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Meme 1 Pa 
[Yn 


My =0 


ri 
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Fig. 5. Torque Memiappi 


bp 


of mechanical losses in the pump with variable capacity qpe.= 


Pn Ap; 


bpp, per one shaft revolution, 


with working liquid reference viscosity v,, as a function of the indicated increase App; of pressure in the pump 
working chambers — graphical interpretation of theoretical models (2) and (7); capacity qp., per one shaft revolution 


(coefficient bp of pump capacity): dpa = 
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Fig. 6. Torque Mpm\4pp;-0,p.v Of mechanical losses in a piston (axial or radial) pump with crankcase filled with liquid and 
with variable capacity dpg= bpp, per one shaft revolution, at the indicated increase App; = 0 of pressure in the pump 
working chambers, as a function of the ratio of viscosity v to reference viscosity v, — v/v,— graphical interpretation of 

theoretical model (4); capacity qr, per one shaft revolution (coefficient bp of pump capacity): qpg= 0 (bp = 0), Are (bp), 

Ire Ie (bp = 1). Torque Mpmappi=0.b,v of mechanical losses in the pump without crankcase filled with liquid is practically 
independent of the liquid viscosity v and is determined at the liquid reference viscosity v, 


AP; 
Mom bp i Npo= cte, Gp, =0 (bp=0),dpgy (bp). Gpg=er (bp=1), Mmins?h,?imax 
App; =0 Ap; 
Mpmlb, = Mpm bp =1 
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Fig. 7. Torque Mrmsppibpv of mechanical losses in a piston (axial or radial) pump with crankcase filled with liquid and 
with variable capacity qpg.= bpqp, per one shaft revolution, as a function of the indicated increase App; of pressure in 
the pump working chambers — graphical interpretation of theoretical model (9); capacity qp., per one shaft revolution 
(coefficient bp of pump capacity): qpg= 0 (bp = 0), deg (be), Ieg= rr (bp = 1); liquid viscosity Vinrin , V, ANA Vinge. Torque 
Mb mdppi.bpv Of mechanical losses in the pump without the crankcase filled with liquid is practically independent of the 
liquid viscosity v and is determined at the liquid reference viscosity v, 


The theoretical and mathematical models describe the 
torque Mpm of mechanical losses in the pump with theoretical 
(constant) capacity qp per one shaft revolution or with 
geometrical (variable) capacity qp,, = bp qp, per one shaft 
revolution: 

pt = QPlAppi = O.ppli = 0.bp = 1vn 1S a theoretical capacity per one 

shaft revolution of the pump with constant capacity per 


one revolution (bp = 1) determined at App, = 0, pp,;= 0 and 
V» Which is equal to the working volume of the working 
chambers created during one shaft revolution, 

drev = bp Gp, is. a geometrical capacity per one shaft revolution 
of the pump with variable capacity per one revolution at 
App; = 0, Pp; = 0 and v,, which is equal to the working 
volume of the working chambers created during one shaft 
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revolution. Capacity qp,, per one shaft revolution changes 
in the 0 < qpa < qr: range and coefficient bp = qp,,/qp, of the 
pump capacity changes in the 0<b,< 1 range. 


The proposed mathematical models describing the torque 
Mbrm Of mechanical losses in the pump, related to theoretical 
models of the torque of mechanical losses, take the form: 

e inapump with theoretical (constant) capacity qp, per one 
shaft revolution (b, = 1) [referring to theoretical model 


(8)]: 
M 


v Ain APpi 
Pm|App;.v ~ k41 Mp, Coy +ky Mp, = 


n n 


Vea AP pj 
= [kai ™ +k4s Aler Mp, = 


n n 


(12) 


v APpi 4 qrip 

=fk awļm ik Pi Ptt n 

I uo 4.2 3 ] oT 
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Mp; pr OP pi 
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M 


Pm|App;. bp=Lv, Mom|app: =0,bp=l, Vn _ 
qdptAPpri 


211 


(14) 


Mom|app =p bp=l, vp Mom|app =0,bp=l, Vn _ 


qdrePn 
2II 


_ Mo mlapp, =p,-bp=L vy, B Momapp, =0,bp=l.v, 2 
Mp, 


e in apump with geometrical (variable) capacity qp,, 
(drev = bp Qr) per one shaft revolution [referring to theoretical 
models (9), (10) and (11)]: 


Manipal (15) 
v App; 
=(k411 +K412 bp)Mp, (m +ky,Mp,bp Poi _ 
n n 


v App; 
=[(Ky 14 +k412 b) oe +ky> bp Pis ]Mp = 


M A i n 
=[(k411 +K412 bp)(—)*™ +k43 bp ae jaeta 
Va Pa 2I 
where: 
Mom|app; =0, bp=0, vy Momlapp =0,bp=0,V_ 
k41 OS 0) 
Mp GptPn 
271 
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Mo ntapy,=0,bp=L, Va —Momlapp;=0. bp=0.Vy _ 


k = 
4.1.2 
Mp i7 
M _ T = M _ o ( ) 
__  Pm|App; =0,bp =L, Vn Pm|App; =0, bp =0.v, 
APPa 
211 
k Sa AM p m|App, bp, Va 
42 7 M 
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Mp 
Commentary: 


— The sum (k,,,+k,,,) of coefficients used in mathematical 
model (15) describing the torque Mpm of mechanical 
losses in the pump with geometrical (variable) capacity 
GPev (Arev = bp qp,) per one shaft revolution is equal to 
coefficient k,, used in the mathematical model (12) 
describing the torque M,,, of mechanical losses in that pump 
working as a pump with theoretical (constant) capacity per 
one shaft revolution: 
kaaa + Kyra = ky. 

— Coefficient k,, used in mathematical model (15) describing 
the torque M,,, of mechanical losses in the pump with 
geometrical (variable) capacity qp,, (Gre, = bp qp,) per one 
shaft revolution is equal to coefficient k,, used in the 
mathematical model (12) describing the torque Mpm of 
mechanical losses in that pump working as a pump with 
theoretical (constant) capacity qp; per one shaft revolution. 


CONCLUSIONS 


1. The theoretical and mathematical models have been 
developed of the torque Mp, of mechanical losses in the 
„working chambers — shaft” assembly of a displacement 
pump with constant qp, (Vp,) and variable qp,., = bp dp; (V pev) 
capacity per one shaft revolution. 

The models describe the relation of the torque Mp,, of 
mechanical losses in the assembly to the torque: 


M, = OPevAP Pi _ bp de: APpi 
Pi 2 OT 


indicated in the pump working chambers and also to the 
working liquid viscosity v at the pump inlet, changing in the 
Vmin <= V < Vmax range. It is assumed that a small change of the 
pump driving motor rotational speed n, (due to the changing 
pump shaft torque M, loading the motor) practically does 
not influence the torque Mpm of losses. 

The indicated torque M,, in the pump working chambers and 
the working liquid viscosity v are parameters independent 
of the torque Mpm of mechanical losses in the „working 
chambers — shaft” assembly. 


The models describe also the relation of torque Mp,, to the 

capacity qp,, per one shaft revolution (coefficient bp = qp../qp; 

of the pump capacity) in a pump with variable capacity per 
onerevolution. 

The assumedchange of qp,, (bp) is in the 0 < qpe < qm 

(0 <b, < 1) range. 

The mathematical models of the torque Mpm of mechanical 

losses are based on defined coefficients k; of energy losses 

relating the torque of mechanical losses to a reference value, 

i.e. to: 

— theoretical torque Mp, of a pump with theoretical 
(constant) capacity qp, per one shaft revolution, 
determined at the increase App; of pressure in the pump 
equal to the nominal pressure p, of system operation 
(App; = Pa), with: 

- known values of the pump capacity coefficient 
bp = dpe/ Apr 

- assumption of practically constant pump speed np 
equal to the speed np, of the unloaded pump shaft 
(np = Npo). 

The mathematical models of the torque Mp of mechanical 

losses in the „working chambers — shaft” assembly should 

correspond with the models of volumetric losses in the 
working chambers and with the models of pressure losses 
in the pump channels. 
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The semi-Markov model of energy state changes 
of the main marine internal combustion engine 
and method for evaluating its operation during 

ship voyage 


Jerzy Girtler, Prof. 
Gdansk University of Technology 


ABSTRACT 


The article presents a method for evaluating the operation of internal combustion engines 
used as the main engines in the propulsion systems of sea-going ships in various operating 
conditions. This method enables calculating the engine operation value based on the theory 
of semi-Markov processes and mathematical statistics. What is noteworthy, the operation 
of the examined internal combustion engines is compared to a physical quantity which 
is expressed in the form of a number with the measure unit called joule-second. a model 
having the form of a semi-Markov process is proposed to describe the energy state changes 


taking place in the main engines during their operation. Also proposed is the use of the point and interval 
estimation at a given confidence level p for determining the value of the energy converted in the engine 
during its operation for the known and unknown standard energy deviation treated as the random variable. 
The semi-Markov model of changes of the energy converted in the engine during its operation is presented 
in a general form. The above model was used for determining the ship 5 main engine operation, which is 
in the examined case a function of the energy converted in particular energy states, the expected value 
of the time of duration of these states and the probability of their existence. These probabilities compose 
the limiting distribution of the semi-Markov process, the values of which are the specified main engine 
energy States. 


Key words: operation; energy; semi-Markov process; piston internal combustion engine 


INTRODUCTION 


References [2, 4, 6, 12, 13] present a proposal to investigate 
the operation of internal combustion engines as the effect of 
the existence of their energy states which makes it possible to 
convert the delivered energy E and transmit it in time t to the 
energy receivers, such as propeller screws, pumps, compressors, 
etc. For the operation understood in the above way a method 
was proposed to valuate the operation of diesel engines used 
in shipbuilding in a general formulation taking into account 
the engine wear [2, 6]. Ref. [2] discusses the operation of the 
abovementioned internal combustion engines used as the main 
engines in the propulsion systems of sea-going ships taking into 
account the fact that correct operation of these engines requires 
securing delivery of sufficient amount of energy to its receivers, 
which are the propeller screws. The method presented in that 
article makes use of the deterministic model of engine operation 
which takes into account ship sailing conditions. 

This article presents a stochastic model in the form of 
a semi-Markov process to model the changes of energy states 
of an arbitrary main engine, and the method for valuating the 
operation of those engines using this model. Like in Ref. [2], 
it was assumed here that the operation (D) can be compared to 
the physical quantity being the product of energy (E) and time 
(t) of its conversion, i.e. D = Ett. 
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The above interpretation of the operation of piston internal 
combustion engines results from the application of the analogy 
method (the intermediate method between induction and 
deduction), which makes it possible to transmit the observations 
from one object of investigation (empirical system) to another. 
Like in Ref. [2], the inspirations for developing the here 
presented method were proposals by P. L. Maupertius and 
W. R. Hamilton to consider the operation of a mechanical 
system as a physical quantity which reflects changes of 
mechanical energy in time. In classical physics an interpretation 
of the operation can be found which considers it the effect of the 
change of energy in time, expressed in the form of the product 
of energy and time, with the joule-second (joulexsecond) as 
the resultant measure unit. 

In this formulation distinction is made between the 
operation [18, 19]: 

e ofthe mechanical system (system of material points) being 
the result of kinetic and potential energy changes, which is 
referred to as the Hamilton operation (D,,) and 

e resulting from the change of only the kinetic energy of 
the mechanical system (body), called the Maupertius 
operation (Dy). 


A similar interpretation of the operation was adopted 
in quantum mechanics with respect to the source of 


electromagnetic radiation [8]. There, the equivalent of the 
operation having a similar sense is the Planck’s constant (h), 
which is also a physical quantity expressed by a number and 
the measure unit [joulexsecond]. 

Achievements in classical physics and quantum mechanics 
in this area were the motivations for the author to introduce the 
operation understood in the above way also to the technique, 
giving it a specific interpretation for particular power generation 
devices, including diesel engines. During the operation of 
these engines the energy conversion also takes place. The 
energy is first converted into heat (Q) and then into work 
(L). In the former case the chemical energy collected in the 
fuel-air mixture is converted into the internal energy of the 
exhaust gas, which bears the name of the thermal energy. This 
process takes place during fuel combustion. In the latter case, 
in turn, the internal energy in the exhaust gas is converted into 
the mechanical (kinetic and potential) energy of the moving 
piston. Obviously, the above energy conversion processes are 
accompanied with losses. 

Such an approach to engine operation seems to be more 
favourable, as comparing energy conversion related advantages 
of particular engines based solely on the analysis and 
evaluation of their power or work does not give comprehensive 
information of the engine capability to perform the task. This 
information can only be obtained from a combined analysis 
of the work and the time of its realisation, represented by 
the quantity D = Lt or D = N-t, as it may happen that the 
engine is capable to convert the required energy to get the 
work L, but the time t which it needs for doing this is longer 
than the required time (for instance, excessively long time of 
duration of transient processes, so-called unsteady states, in 
the engine before it is fully loaded). Or, it may happen that 
the work L = Nt needed for performing the assumed task will 
not be done by the engine in a given time t (for instance, due 
to excessively small power output N of the engine resulting 
from its wear). Hence the conclusion that evaluating the 
correctness of engine operation based only on its power or 
work, without simultaneously taking into account the time of 
engine operation, is insufficient for full assessment of its energy 
conversion related qualities. Obviously, the engine power is an 
important parameter, as it contains the information how quickly 
the work can be done. But the operation in the above presented 
formulation is also important as it contains the information how 
long the engine can work. Since obtaining the required work 
(L) via conversion of the energy (E) is accompanied by energy 
losses in the form of heat (Q), therefore it is more reasonable 
to analyse energy changes, and not only engine work during 
engine operation. 

The advantage of the interpretation of the diesel engine 
operation proposed by the author is that the earlier descriptive 
evaluation of engine operation, for instance: the operation is 
good, bad, etc., can be replaced by the assessment based on 
comparing the engine operation with the reference standard 
operation using numbers and the measure unit, which is the 
joule-second. 

The sense of the above interpretation of the operation of 
an engine, and in general any power generation device, can 
also be justified by the observation that changes of motion of 
a body depend on the force (F) acting on this body and the time 
(t) during which this action takes place. Body’s capability to 
move is expressed by the product of force and time (F-t) bearing 
the name of the impulse of force [9]. The measure unit of the 
impulse of force is the newton-second. By analogy, we can 
conclude that the operation (D) of the diesel engine depends 
on work (L) done by this engine and time (t) of its operation, 
Le. D=Ltt. 


CONDITIONS OF MAIN ENGINE 
OPERATION AND THEIR STOCHASTIC 
NATURE 


The conditions of main engine operation depend on the 
conditions in which it works, which are determined by external 
conditions of sailing of the sea-going ship and tasks undertaken 
by ship users (the crew). These conditions and tasks are the 
reason why different amounts of energy are converted in 
different times in the working spaces of these engines. In each 
case, however, the operation of the main engine is defined by 
the area of its performance [2, 11, 16, 17]. These performance 
areas are determined by engine speed characteristics, including 
external and control characteristics. If these characteristics 
are projected against the propeller characteristics, which also 
belong to the speed characteristics, then the ranges of operation 
of these engines are defined [4, 5]. 

It was shown in Ref [2] that the points of cooperation 
for a given propeller characteristic are created depending 
on the injection pump setting. Therefore it may happen that 
along a given propeller characteristic the main engine can be 
loaded with power according to the following external power 
characteristics [11, 16, 17]: 

e external characteristic of partial power N.. (h, = idem, 
c=1,2,...,N..< Neu) 
e external characteristic of continuous operating power N „y 

(h, = idem), 

e external characteristic of nominal power N.a 
e external characteristic of maximal power N 
= idem). 


(h, = idem), 
(Hie = 


emax max 


A sample realisation of power changes of the main engine 
during its operation is shown in Fig. 1 [2]. The effective 
power N, generated by the main engine characterises the 
engine operation in the aspect of the rate of energy conversion 
into work, taking into account various types of losses, in 
particular thermal loss. However, engine operation consisting 
in energy conversion into work is not possible if this energy 
is not earlier converted into heat in the engine working spaces 
[2, 15]. Therefore when analysing the main engine operation we 
should take into account the energy delivered to the engine in 
the fuel-air mixture, which is first converted into heat (Q) and 
then into work (L), rather then the pure engine power output. 


tc tcp tb toe te 


ta tag tB 


Fig. 1. Sample realisation of the process of main engine load changes 
during engine operation: N, — effective power, t — operating time [2] 


The analysis should focus on the process of changes of 
the energy delivered to the engine. Particular states of this 
process should be the energy states which release the energy 
Ej = 1, ..., 4), and secure the generation of particular powers 
N,G= 1, ...,4), which are essential for correct engine operation 
and secure realisation of the operating task, (Fig. 1): 

Ni = Nao N, =N 


ec? etr? N; = Nens N, = Naia (1) 
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The analysis can be performed using a deterministic or 
stochastic model for describing real main engine energy change 
processes. 

The deterministic model of energy state changes taking 
place in the main engine during its operation was presented and 
discussed in Ref. [2]. This article presents the stochastic model 
of these changes having the form of a semi-Markov process. 


MODEL OF ENGINE OPERATION IN THE 
FORM OF A SEMI-MARKOV PROCESS 


Developing the model of internal combustion engine 
operation in the form of the semi-Markov process requires 
the application of the theory of semi-Markov processes. The 
publications on the semi-Markov processes give definitions 
of this process which differ by the range of generality and 
precision [3, 7]. For the purpose of modelling changes of the 
internal combustion engine energy states the semi-Markov 
process (family of random variables) {X(t): t > 0} at T = [0, 
+0), can be defined using a so-called homogeneous markovian 
process of recovery [7]. 

The semi-Markov model of an arbitrary real process can be 
only constructed when the states of this process can be defined 
in such a way that the time of duration of the state existing in 
time t, and the state which is to be obtained in time T, do not 
depend stochastically on the states which earlier took place, 
nor the time intervals of those states. 

Constructing the semi-Markov model {X(t): t > 0} of the 
real process of energy state changes taking place in the diesel 
engine operation phase is possible, because [3]: 

1) the Markov condition is fulfilled which says that the future 
evolution of an arbitrary process of the energy state changes 
taking place during engine operation depends solely on 
the engine state at a given time instant and not on the past 
engine operation, consequently the future state of this 
process and the time of its duration depend on the present 
and not on the past, 

2) random variables T, which represent the time of duration 
of the state e; irrespective of the fact which state will be the 
next, and T;,, which represent the time of duration of the state 
„e on condition that the next state of this process is the 
state ,,¢;”, have distributions different than the exponential 
distribution. 


When analysing changes of the main engine energy states 
e; resulting from converting particular energies E,(i= 1, ..., 4) 
in the interpretation expressed by formula (1) for engine 
operation which enables the ship to perform the transporting 
task in a relatively long time, theoretically tending to infinity 
(t — œ), we can omit time intervals connected with changes 
of energy states e, taking place when the energy converted in 
the engine working spaces increases from E, to E,,, to generate 
engine power output corresponding to those energies (1) in 
given time intervals (Fig. 1). Consecutive power outputs 
N, G = 1, ..., 4) can be considered the effects of the engine 
energy states e; € Ej = 1, ..., 4) of the stochastic process 
{X(t): t> 0}, each of which takes a constant values in the time 
interval [Tns Tm). 

That means that the semi-Markov process {X(t): t > 0} 
is the process having the realisations constant in intervals 
and continuous on the right. Therefore the process is discrete 
in states and continuous in time. a sample fragment of 
the realisation x(t) of this process is shown in Fig. 2. This 
realisation reflects the situation when, following the elapsed 
time t, successive state changes of this process take place in 
successive time intervals [T,,, Tm). 
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Fig. 2. Fragment of a sample realisation x(t) of the semi-Markov process 
(X(t): t= 0} 


In case of marine main engines, in an arbitrary time of 
their operation the process {X(t): t > 0} can take one of the 
energy states existing during engine operation and belonging 
to the set: 


B= {ees €z, €43} (2) 


with elements having the following interpretation (Fig. 1): 
engine energy state resulting from generation of energy 
E, corresponding to engine load with partial power 
N, = Neo 

engine energy state resulting from generation of energy 
E, corresponding to engine load with (continuous) 
operating power N, = Naw 

engine energy state resulting from generation of energy 
E, corresponding to engine load with nominal power 
N; = Nan 

engine energy state resulting from generation of energy 
E, corresponding to engine load withy maximal power 
N, = Nois 

Making distinction between those energy states differing 
by engine load is important because these states originate from 
the principles of diesel engine operation and are connected with 
the need for precise positioning of the fuel lever corresponding 
to particular loads [11, 15, 16, 17]. 

According to the presented model of energy state changes in 
the form of the process {X(t): t > 0}, an arbitrary main engine 
which is in state e, during its operation can reach, due to load 
increase, the state e, with the probability p,, after time T,,, and 
then reach the state e, with the probability p,, after time T». 
The state e, of the process {X(t): t= 0} will take place when an 
arbitrary main engine starts converting energy E,. Ina similar 
way, from the state e, the process can reach the state e, with 
the probability p,, after time Tą, or the state e,, which can take 
place with probability p,, after time T,,, etc. 

A characteristic feature of the presented model of main 
engine energy state changes {X(t): t > 0} is that the change of 
the state e; into the state e; only depends on the properties of 
the state e;, and not on the previous states. This feature justifies 
the opinion that it is the semi-Markov process. 

In this process the set of states E, = {e,, €5, €, €4} is finite. 
The elements composing this set are the values of the semi- 
Markov process X(t): t > 0}. Changes of states in this process 
take place at times Tọ = 0, T,, T>, ... and after the lapse of time 
intervals of their duration (Fig. 2) being random variables with 
different distributions. The state changes in this process take 
place with the following probability: 

P{X (tau) = ejs Thsi Tn < (3) 


<1 X(t.) Hey, RC eave X Ta — Ta» Tja to} = 
= PAC) = €j Th41 — Tn < tX) =e;} 


where: 
e, & € E,i,j=1, 2, 3, 4; e Fe, 


n+l 


As a result, when the state ei = 0, 1, ..., 4) of the process 
{X(t): t > 0} is known at time qt, the time of its duration and 
the state e; reached at time 1,,, do not depend stochastically on 
the process states taking place at times 0 = Tp, Ti, ..., Tiy , Nor 
the time intervals of their duration. 

The stochastic process {X(t): t >0} is the process whose states 
are identical (constant) in intervals and the realisations of these 
states are continuous on the right. The lengths of the intervals [t,, 
T1), [Ti T2), [T2 T3), -. -> [Tns Tres)» -< in Which the process {X(t): t 
>0} takes constant (the same) values are random variables with 
positive distributions. In case of the loads of the abovementioned 
engines we can assume [3] that the duration time of an arbitrary 
state e; E E. (6), which was reached at time t, and the state 
reached at time t,,, do not depend stochastically on the states 
which earlier took place, nor the time intervals of their duration. 
Consequently, the process {X(t): t>0} is the model of randomly 
changing main engine loads, with the set of states E, = {e}, €», €3, 
e,} and the graph of state changes shown in Fig. 3. 

Pi2 P23 P34 


Fig. 3. Graph of state changes in the process {X(t): t 20} 


The specified states e; € E, (i = 1, 2, 3, 4) of the process 
{X(t): t>0} of an arbitrary main engine can be recognised using 
relevant diagnostic systems (SD), the applicability of which 
depends on the quality of the adopted diagnosing system (SDG) 
and its capability to recognise the states of the abovementioned 
engines as the diagnosed systems (SDN). 

During the operation of an arbitrary diesel engine working 
as the main engine, changes of the states belonging to the 
set E, = {e; 1 = 1, 2, 3, 4} can be interpreted as the process 
{X(t): t = 0} with the realisations constant (identical) in 
consecutive time intervals and continuous on the right [3, 7]. 
The lengths of the intervals in which the process {X(t): t> 0} 
takes the constant (identical) values are the random variables 
T; representing the time of duration of the state e; € E of this 
process, on condition that the next state is e; € E. where 
i,j=1,...,4andi#j. These variables are independent random 
variables with finite expected values E(T;,) and have positively 
concentrated distributions. Moreover, this process reveals 
a property consisting in the fact that the duration time of the 
state e; which took place at time qt, and the state which started at 
time T, do not depend stochastically on the states which took 
place earlier, nor the time intervals of their duration. Therefore 
we can assume that the future states (situations) depend solely 
on the presently existing situation. That means that the process 
{X(t): t > 0} is the semi-Markov process with the state change 
graph shown in Fig. 3. Defining this process requires defining 
its initial distribution P, and the functional matrix Q(t). 

The initial distribution of the process {X(t): t > 0} is the 
following: 


1 for i=0 
P; = P{X(0) = e;} = (4) 
0 for i=1,2,...,3 


According to the state change graph shown in Fig. 3, the 
functional matrix has the following form: 


0 Qt 0 0 
0 Qa t) 0 Q33 (t) 
0 0 Q39(t) 0 


The elements of the matrix (5) are non-decreasing functions 
of the variable t which represent the probabilities of switching 
the process {X(t): t > 0} from the state c; to the state e; 
(e, § € Es i,j = 1, 2, 3, 4; 1 £j) in time not longer than t, and 
are described in the following way [3, 7]: 

Q(t) =P Xa = 6 
= ej Tui 1, <t | X(t) = e} = piF;,(t) 
where: 
e, & € E.G, j= 1,2, 3,4;i +j), 
and: p; is the probability of switching in one step of the 
homogeneous Markov chain inserted in the process {X(t): t> 0}. 
F,(t) is the distribution function of the random variable T; 


representing the duration time of the state c; of the process 
{ X(t): t= 0} on condition that the next state will be e.. 


The probability p, is interpreted in the following way: 
Pi = PAX) = 6) X,) =e} = imQ 
In this situation solving the problem consist in finding the 


limiting distribution of the process {X(t): t > 0} having the 
following interpretation: 


P= limP{X(t)=e,}, j=1.4 
j t>0 
This distribution can be calculated using the formula [7]: 


nm E(T;) 


P = ,j=b234 @ 
> TE(T, ) 
k=0 
where: 
x, = lim tY PXC) =e,|X(0) =e} 
noon yy 


[x,;j= 1, 2, 3,4] — the stationary distribution of the Markov 
chain {X(t,): n € N} inserted in the 
process {X(t): t> 0}. 


The matrix (5) is the stochastic matrix, therefore this 
distribution is fulfilled by the following equation system 
(19) [7]: 


0 1 0 0 
Pa 9 P3 = 
[n T, T3, T4] = [n]; T3, 03,74] 
0 Pa 0 Pq (9) 
0 0 1 0 


ntm tmr = 1 


After solving the equation system (9) we obtain the 
following relations according to the formula (8): 


P= PoiPsE (Th) . P = Par ad, 


2 
M M (10) 
pre PEC), P = P23P34E (Ty) 
3 M 2 *4 M 
where: 
M = pop3E(T,) + p3,E(T,) + Pp3E(T;) + PosP34E(T,); 
Pp; | — probability of switching the process {X(t): t> 0} from 


state e; to state e; (e, e, € E3 i,j = 1, 2, 3, 4; 1 #j); 

E(T;) — expected value of the random variable T,(j = 1, 2, 3, 4) 
representing the duration time of the state e; € E, 
(j= 1, 2, 3, 4) of the process {X(t): t= 0} irrespective 
of the next state to which this process switches. 
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The expected values E(T;) are related with the expected 
values E(T;,) and the probability p, in the following way: 


E(T;) =E(T,) =X pE Lj LA i435 ay 
j 


Therefore, according to the adopted interpretation of the 
internal combustion engine operation [2, 10], the main engine 
operation can be defined in the following general way, see (Fig. 4) 
and formula (7): 


Dhists)= PE,E(T,) + P,E,E(T,) + PEE ;)+ 
+ P,E,E(I,) + P,E,E(L,)+P,E,ECL)+ (12) 
+ PE E(T,) + P,E,E(T) 


The main engine operation defined by formula (12) takes 
place when we can assume that the changes of the propeller 
characteristic, resulting from changes of sailing conditions, do 
not affect considerably the changes of energy conversion in 
the engine working spaces, i.e. we can assume that the energy 
E = idem. 

When changes of the energy E are larger and cannot be 
omitted, the formula (12) is to be modified by introducing the 
average value of this energy, which, as the random variable, 
is the statistics E. Then the operation will be defined in the 
following way: 


Dkr. t )= PEET) + P,E,E(T,) + P,E,E(T,) + 
+ P,E,E(T,) + P,E,E(T) + P|E;E(T;) + 
+ PE E(T,) + P,E,E(1) 


A more accurate definition of the operation requires 
introducing probabilities P,( = 1, ..., 4) given by formula (20) 
to equations (12) and (13). 

Particular probabilities P; = 1, ..., 4) have the following 
interpretations: 


P, = limP{ X(t) = e,}, P, = limP{X(1) = e3}, 


(13) 


r= limP {X = e}, P, = limP {X(t) = e,} 


The probability P, can be interpreted as the probability of 
engine load according to the partial power characteristic, and 
the three remaining probabilities P,, P,and P,can have similar 
interpretations referring to the operating power, nominal power, 
and maximal power characteristics, respectively. 

Obtaining the (obviously approximate) values of the 
probabilities P; = 1, 2, 3, 4) requires valuating p, and E(T)). 

Valuating the probabilities p; and the expected values E(T;) 
is possible when we have the realisation x(t) of the process 
{X(t): t= 0} ina relatively long time interval of investigations, 
i.e. for t € [0, t,], where t, >> 0. Then we can obtain numbers 
n,(i,j = 1, 2, 3, 4; 14j), which represent the number of changes 
from the state e; to the state e; in the relatively long time. 

The estimator of the highest likelihood of the probability 
of change p; is the statistics: [7] 


Ê. =-=}, i#j; ij=1, 2, 3, 4 


(14) 


the value: 


2 Ni 
j 
of which is the assessment of the unknown probability p, of 


switching the process {X(t): t > 0} from the state e; to the 
State e. 
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From the abovementioned realisation x(t) we can also 
obtain realisations t,(™, m = 1, 2, ..., n; of the random variables 
T;. The application of point estimation provides opportunities 
for easy assessment of E(T;) as the arithmetic mean of the 
realisation t,(™. 

The expected value of energy in formula (12) is the observed 
value of the statistics E,(k = 1, 2, 3, 4). Obviously, this statistics 
is the random variable [1, 7] having the following general 


formula: 
— 1 n 
N iz 
where: 
E; — random variables having the same (arbitrary) 


distributions, with the same expected value E(E,;) = 
=m), and variance: 


2 2 
D (Eki) = Oki #0 
This way the random variable E, (statistics) is the arithmetic 
mean of n independent random variables E, having the same 


distributions. The expected value and the variance of the 
random variable E, are given by the relations [1, 14]: 


E(E,) = E(E,) = Miks 


r a 82 (16) 
D? (Em) = >D? (Ey) =% 
n n 

It results from the Lindeberg-Levy theorem [1] that the 
random variable E, has the asymptotically normal distribution 
N(m,,, 6,/VN) irrespective of the nature of the random variable 
E,. That means that the arithmetic mean of n independent 
random variables E, having an arbitrary but the same 
distribution and the same expected value, E(E,;) = m,,, and the 
variance D2(E,;) = of, has the asymptotically normal distribution 
N(m,,, o Vn). 

Applying the proposed methodology to the analysis of 
energy changes during the main engine operation seems 
attractive due to the fact that the convergence of the distribution 
of the statistics E, to the normal distribution N(m,,, o,/VN) is 
very fast and we can use it for all n > 4, i.e. always in practice 
[1]. 

When the value o, is known, taking into account the 
distribution N(m,,, 6,/VN) of the statistics E, we can calculate 
the confidence interval for the unknown expected value m,, = 
E(E,) using the formula [1, 10]: 


Pla, =A + <E(E,) SE +y, S| =$ 07) 


where 

B -— confidence level, 

y, — standardised variable of the normal distribution, 
which corresponds to the confidence level B = 1 — a 
(a - significance level), 

Ek — value of the statistics Eo which is the arithmetic mean 


of the performed measurements. 


In experimental practice, as a rule, the value o, is not known, 
but it can be evaluated using the formula [1, 14]: 


S= YG -& (18) 
i=l 


Since the statistics E, has always the asymptotically normal 
distribution N(m,,, 6,/VN), and its convergence to the normal 
distribution N(m,, ©) is very fast, we can assume in practice that 


the investigated energy E, has normal distribution N(m,,, 6,). 
It is noteworthy, however, that adopting the above assumption 
is not a limitation in experimental practice, as the statistics E, 
has always statistically normal distribution N(m,,, 6,/Vn). 
Moreover, the convergence of this distribution to the normal 
distribution is very fast [3]. Therefore the random variable: 


Ey - Hes) (a4 
Sk 
has the t-Student distribution with k = n — 1 degrees of 
freedom [1, 10]. That means that the confidence interval for 
the unknown expected value E(E,) of the random variable E, 
can be calculated using the formula [1, 10]: 


S S 
Piety n < FEY) < Ektta p= = 8 09) 
g a, | (Ep k" ‘a, 5 B 


In practice, the formula (19) is more applicable for assessing 
the energy changes during engine operation in real conditions 
of ship sailing than the commonly used formula (15), as it not 
only reflects the randomness of the results obtained in operating 
conditions, but also determines the error in estimating the 
expected value of the energy converted in the engine in its 
particular states. This error, which in the present case is equal 
to 1 — B, is the not covering of the expected value E(E,) by the 
confidence interval given by the formula (19). 

The application of the estimation in intervals makes it 
possible to determine the expected value E(E,) in the form of 
an interval with random limits, which means that the following 
inequality is to be fulfilled with probability B: 


E,<E(E,) <E, (20) 


It results from the inequality (20) that the main engine 
operation which can be considered safe can be assessed in the 
following way, according to formula (13): 


Dit, t9) = PE, gE(T,) + P,E24E(T,) + P3E3gE(T3) + 
+ P E44E(T3) + P,E,gE(T,) + P3E3gEC3) + 
+ PE 4E(T,) + P,E.gE() (21) 
On the other hand, the optimistic engine operation can be 


assessed, according to formula (13), as: 


D[t,,t9)= P/E; ,E(T,) + P E2E(T3) + P3E3.E(T3 ) + 
+ PyE ygE(Ty) + PyExgE(Tp) + PsEygE(T3) + 
+ PE) E(T,) + P,E,,E(T)) (22) 


If necessary, processes of this type with a larger number of 
states e, 1= 1, 2,3, ...K can be investigated. Taking theoretically 
into account these needs, we can most generally discuss the 
process {X(t): t > 0} with the set of states E,, = {€), &, €3, ..., 
e,}. This process will be fully defined if, like in previous cases, 
a formula is given for the functional matrix: 


Q; -_ [Q; (Oli = 123; s K 
and the initial distribution: 
P.=P{X(0) =e} i=1,2,3,...,K 


For this generalised case the matrix of the process 
{X(t): t > 0} is the following: 


0 QW 0 an 0 0 0 
Q(t) 0 Q(t)... 0 0 0 


(23) 


Q(t)= (24) 


0 0 O Qkit) 0 Qkit) 
0 .. 0 Qk) 0 


and its initial distribution has the form: 
1 for i=1 
P= P{X(0)=e)} = 


(25) 
0 for i= 2,3,....K 
Consequently, P{X(0) =e,} = 1. 


Like in previous analyses, we have to assess the limiting 
distribution of the process [7]: 


nt E(T; 
p= (a 2. seg K 


J K 
SmE) 
k=1 


Using the same procedure as in previous cases we arrive 
at the relation: 


(26) 


Il a E Pk-ik-2) E(T;) 


oa U j=2,3,..., K 
l Sy A-Pkaik) 
ET) +| X | [—©* EC) (27) 
j=2 k=2 kk-1 


It results from formula (27) that, for instance, the probability 
of engine load resulting from the appearance of the state e, can 
be calculated as [7]: 


E(T, 
P, —r r (28) 
E(T,) + TL = Pk-ik-2) E(T)) 
j=2k=2 kk-1 


The presented considerations reveal that after adopting the 
interpretation of the operation in the here proposed version and 
making use of the theory of semi-Markov processes we can 
determine the operation of marine main engines in probabilistic 
formulation, which is more adequate than the deterministic 
formulation, as it reflects random conditions of operation of 
engines of this type. 


REMARKS AND CONCLUSIONS 


- The article presents the probabilistic method for evaluating 
the operation of an internal combustion engine used as 
the main engine on a ship. This method can be applied 
to an arbitrary piston internal combustion engine of both 
spark-ignition and compression-ignition type. The method 
takes it into account stochastic nature of energy conversion 
processes observed during engine operation. It was proved 
that a very useful model for investigating those processes 
is the model having the form ofa stochastic process which 
is discrete in states and continuous in time. The analyses 
presented in the article reveal that this model can be the 
semi-Markov process which has the abovementioned 
features. 

- The operation of the internal combustion engine was 
interpreted as generation of the required energy in a given 
time and its delivery (transmission) to the receiver. That 
means that when analysing the operation of these engines 
in given time and taking into account both the energy 
generated by these engines and time of its generation, we 
can compare (in a valuating formulation) this operation to 
a physical quantity which has a numerical value and the 
measure unit called joule-second [joulexsecond]. 

When analysing energy related qualities of internal 
combustion engines we should analyse their whole 
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operation, and not only their work, as the engine operation 5. Girtler J.: Energy-based aspect of operation of diesel engine. 

includes not only the energy conversion into work, but also COMBUSTION ENGINES No 2/2009 (137). 

into heat. 6. Girtler J.: Conception of valuation of combustion engine 
operation. Journal of KONES. Powertrain and Transport. 
Editorial Office Institute of Aeronautics BK, Warsaw 2008. 

7. Grabski F.: Theory of semi-Markov processes of operation of 


- The process of changes of energy states during the operation 
of an arbitrary diesel engine is a random process with 


continuous, positive and limited realisations. technical objects (in Polish). ZN AMW nr 75A, Gdynia 1982. 
- The models of changes of energy states may refer to the g, Gribbin J.: In Search of Schrödinger’s Cat Quantum Physics 
operation of an arbitrary engine, and include various Reality. (Polish issue). Zysk i S-ka Wydawnictwo s.c. Poznań 
numbers and interpretations of energy states. In each case, 1 997. 
however, a model can be constructed in the form ofa semi- 9. Helwitt P.G.: Physics around us (in Polish). PWN, Warszawa 
Markov process, continuous in time and discrete in states, 2001. 
i.e. having a limited number of states. 10.Pawłowski Z.: Mathematical statistics (in Polish). PWN, 
- The developed semi-Markov model of the process of Warszawa 1980. 


11.Piotrowski I., Witkowski K.: Operation of marine internal 
combustion engines (in Polish). AM, Gdynia 2002. 
12.Rostanowski J.: Identification of ships propulsion engine 


energy state changes during the main engine operation is 
the process with limited number of states and continuous 


ee : Tani operation by means of dimensional analysis. Journal of POLISH 
Another approach to investigating the energy converted CIMAC, Vol 4, No 1, 2009. 

in the main engine and transmitted to the propeller screw 13 Rudnicki J.: On making into account value of operational 

was also proposed. This approach consists in considering applied to ship main propulsion engine as an example. Journal 
this energy a random variable. The applicability of the of POLISH CIMAC, Vol 4, No 1, 2009. 

Lindeberg-Levy theorem was proved which says that the 14.Volk W.: Applied statistics for engineers (in Polish). WNT, 
statistics created based on the results of the investigations Warszawa 1965. 


of a given energy is the random variable having the 15.Wajand J.A.: Compression-ignition engine (in Polish). WNT, 
asymptotically normal distribution, irrespective of the Se oR À ce 
nature of the random variable representing this energy. This 16.Wiodarski J.K.: Operating states of marine internal combustion 


À fh ae co engines (in Polish). WSM, Gdynia 1998. 
enables using the estimation in intervals to evaluate the 17.Wojnowski W.: Marine internal-combustion power plants, Pt I. 


unknown expected value of the random variable, which can (in Polish). AMW, Gdynia 1999. 
be the energy converted in the main engine and transmitted 18.Encyclopaedia of contemporary physics (in Polish). Collective 
to the propeller screw during engine operation. work. Editor: Redakcja Nauk Matematyczno-Fizycznych i 
- The application of the estimation in intervals to the Techniki Zespołu Encyklopedii i Slownik6w PWN. PWN, 
evaluation of engine operation has made it possible to assess Warsaw 1983. 
its value in an optimistic version (22) and a safe version 19.Scientific and technical lexicon with supplement (in Polish). 
(21). Collective work. Editor: Zespół redaktorów Działu Stownictwa 
Technicznego WNT. WNT, Warsaw 1989. 
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ABSTRACT 


The article compares combined systems in naval applications. The object of the analysis is the combined 

gas turbine/steam turbine system which is compared to the combined marine low-speed Diesel engine/steam 

turbine system. The comparison refers to the additional power and efficiency increase resulting from the 

use of the heat in the exhaust gas leaving the piston engine or the gas turbine. In the analysis a number 

of types of gas turbines with different exhaust gas temperatures and two large-power low-speed piston 

engines have been taken into account. The comparison bases on the assumption about comparable power 
ranges of the main engine. 


Key words: marine power plants, combined systems, piston internal combustion engine, gas turbine, steam turbine 


INTRODUCTION 


In recent years combined systems consisting of a gas turbine 
and a steam turbine have started to be used as marine propulsion 
systems. In inland applications the efficiency of these systems 
can exceed 60%. In naval applications such a system was 
used in Millenium, a passenger liner. But it is, so far, a unique 
case of propulsion system related with the highest efficiency. 
However, this system requires a more expensive fuel, which 
is the marine Diesel oil. 

Another possible combined system can consist of a piston 
engine which cooperates with the steam turbine utilising the 
heat in the exhaust gas leaving the piston engine. The leading 
engine in this is system is always the piston internal combustion 
engine. The objects used as the main engines will be large 
low-speed piston engines which burn heavy fuel. At present, 
the efficiency of these engines reaches 45-50%. In case of such 
large main engine power ranges, the exhaust gas leaving the 
piston engine contains huge amount of heat which can be the 
object of further utilisation. 

Another possible solution of a combined system which 
can be applied in shipbuilding consists of a gas turbine, and 
a steam turbine which utilises the heat leaving the gas turbine. 
This heat can be used for production of steam which, in turn, 
can be used in the steam turbine cycle. 


CONCEPT OF A COMBINED SYSTEM 


Combined propulsion systems in naval applications are first 
of all used in fast special ships and in the Navy as the systems 


being the combination of a Diesel engine and a gas turbine 
(CODAG, CODOG) or gas turbines (COGOG, COGAG). The 
propulsion applied in the passenger liner Millenium makes 
use of a COGES system which increases the efficiency and 
operating abilities of the ship by combining the operation of 
a gas turbine with a steam turbine which drives the electric 
generator, while the propeller screws are driven by electric 
motors. In this system the steam turbine cycle is fed with the 
steam produced in the waste heat boiler supplied with the 
exhaust gas from the gas turbines. 

Fig. 1 shows sample efficiency curves of the combined gas 
turbine/steam turbine systems as functions of power plant load, 
as compared to the gas turbine working in a simple open cycle 
and the marine low-speed Diesel engine. The analysis of the 
efficiency changes reveals that the combined system composed 
of a gas turbine and a steam turbine has the highest efficiency, 
which can reach up to 60% for maximal loads. Gas turbines 
working in a simple open cycle have the lowest efficiency, 
of an order of 33+35% on average, and about 40% as the 
maximum [7]. Low-speed Diesel engines reach efficiency of 
an order of 47+50% [5, 6]. At the same time their efficiency 
characteristics as functions of load are flat, which is of highest 
importance in marine propulsion systems working in heavily 
changing load conditions. For instance, the relative efficiency 
drop An/neris equal to 15+20% when the load changes from 
100% down to 50% in the combined gas turbine/steam turbine 
systems and gas turbines working in a simple cycle, while for 
a low-speed Diesel engine this efficiency drop is equal to 1+2%. 
This property of the Diesel engine and utilisation of additional 
heat contained in its exhaust gas makes the application of this 
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system very useful for marine propulsion systems working in 
extremely changing load conditions. 

The temperatures of the exhaust gas from the gas turbines 
are equal to 450+600°C on average, while for the low-speed 
Diesel engine these temperatures are of an order of 220+300°C 
[5, 6], Fig. 2. But in the gas turbines this temperature remarkably 
decreases with decreasing load, while in the Diesel engine its 
changes are not high, nor monotonic: initially the temperature 
decreases and then starts increasing with the decreasing load. 

This property of the Diesel engine in partial loads makes it 
possible to keep the temperature of the live steam at a constant 
level within a wide range of steam turbine cycle loads. 
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Fig. 2. Temperatures of the exhaust gas from the marine Diesel engine 
and the simple gas turbine as functions of power plant load 


Power evaluation 
of the combined propulsion system 


The power of the combined propulsion system is 
calculated by adding up particular power outputs of the system 
components, i.e. the main engine and the steam turbine: 


N combi = Ne + Nsr (1) 
while the efficiency of the combined system is: 
n = Neombi =np: 1+ Nor (2) 
combi my: Wu D Ne 
where: 
Nc Mp — are the efficiency and mass flow rate of the main 


engine. 
It results from the relation (2) that each additional power 


in the propulsion system increases the efficiency of the system 
and, consequently, decreases the fuel consumption. The 
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higher the additional power obtained from utilisation of the 
heat contained in the main engine exhaust gas, the lower the 
fuel consumption. In these circumstances we should tend to 
reach maximal power also from the steam turbine, in which 
the power output largely depends on proper selection of live 
steam parameters and condenser parameters. 


Steam turbine cycle 


The combined cycle makes use of the waste heat in the 
exhaust gas from the Diesel engine or the gas turbine. Adding 
a steam cycle to the combined Diesel engine or gas turbine cycle 
makes it possible to increase the power of the combined system, 
i.e. increase the system efficiency according to formula (2). 

In case of small powers and low live steam temperatures in 
combined systems, single pressure systems are used [1, 2, 3], 
Fig. 3. Such a system consists of a single pressure waste heat 
boiler), a condensing steam turbine, a water cooled condenser, 
and a single stage feed water preheater in the deaerator. 


Fig. 3. Flow diagram of a single pressure system 


The application of the single pressure system does not secure 
optimal utilisation of the energy in the exhaust gas when its 
temperature is high. The system which is most frequently used 
in these cases includes an additional low pressure evaporator 
[2, 3], Fig. 4, which not only increases the utilisation of the 
waste heat in the exhaust gas, but also provides opportunities 
for better thermodynamic use of the low-pressure steam. 


Fig. 4. Flow diagram of a two-pressure system 


Limits of steam cycle parameters 


Limits of steam cycle parameters result from strength, 
technical, and durability conditions for particular system 
elements, but also from constructional and economic 
restrictions. 

The temperature difference between the exhaust gas and 
the live steam was assumed equal to At = 10°C in utilisation 
boilers used in shipbuilding, following the literature [3, 4]. The 
“pitch point” value recommended by MAN B&W [4] for marine 
boilers was assumed equal to ôt = 10°C. The limit of the steam 
dryness factor x behind the steam turbine was assumed equal 
tO Ximi = 0.86-0.88. For marine condensers cooled with salt 
water, the company MAN B&W [4] recommends assuming the 
condenser pressure equal to pg = 0.065 bar. The temperature of 
the water which feeds the boiler is of high importance for the 
lifetime of the feed water heater in the boiler. MAN B&W [4] 
recommends assuming that the feed water temperature is not 
lower than 120°C, when the content of sulphur exceeds 2%. 
The caloric value of the fuel was assumed constant and equal 
to Wu = 42700 kJ/kg in the calculations. 


Steam cycle optimisation 


The combined systems, in this case the steam cycles, should 
make maximal use of the heat contained in the exhaust gas 
from the main engine. Therefore the optimisation is reduced 
to calculating steam cycle parameters for which the steam 
turbine reaches the maximal power output. The area within 
which optimal steam cycle parameters are searched for is to 
be narrowed to the area in which all above listed limits and 
restrictions imposed on the steam cycle are fulfilled. 


In case of combined systems with gas turbines, a two- 
pressure steam turbine cycle is used. The calculations for the 
combined system with Diesel engines were performed for both 
the single- and the two- pressure variant. The steam turbine 
cycles in the combined systems were calculated based on the 
same schemes, and the same above listed assumptions and 
parameters. 


COMBINED GAS TURBINE/STEAM 
TURBINE SYSTEMS 


In the calculations of the combined systems with gas 
turbines a number of types of gas turbines produced by 
Alstom and GE were taken into account. The analysed turbines 
represented different power outputs and different exhaust 
gas temperatures. The total power output of the gas turbines 
in the combined system was kept close to each other in all 
variants. The gas turbine parameters which were assumed in 
the calculations are collected in Table 1. 


Tab. 1. Parameters of gas turbines used in the combined 
gas turbine/steam turbine cycle [7] 


Type 
of turbine 


CYCLON 
GTX100 
GT10B 
LM2500 
GT35 


Tab. 2. Combined gas turbine/steam turbine cycle with two-pressure boiler 


Type CYCLON 


GTX100 | GT10B | LM2500 


Number of engines 


g/kWh 


kg/kWh 


% 


% 


kg/kWh 
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In this case the steam turbine system was analysed in the 
cycle with the two-pressure boiler, see Fig. 4, for the feed 
water temperature tpw = 120°C. The results of the combined 
cycle calculations are given in Table 2. The power output of 
the steam turbine with respect to that of the gas turbine changes 
between 35% and 49%. The power output of the steam turbine 
increases mainly due to the increase of the gas turbine exhaust 
gas temperature, but also depends on the ratio of the mass 
flow rate of the gas turbine exhaust gas to the turbine power 
output. The live steam pressure depends on the gas turbine 
exhaust gas temperature and increasers with the increase of 
this temperature. The pressure p, of the low-pressure stage also 
increases with the increasing exhaust gas temperature (t,), but 
this increase is not large. 

The mass flow rates m, of the heating steam taken from 
the turbine extraction point for heating the condensate in the 
deaerator are relatively comparable. The mass flow rate m, of 
the low-pressure steam which is passed from the boiler to the 
steam turbine decreases with the increasing temperature of the 
gas turbine exhaust gas. 

The combined cycle applied to a system with the gas 
turbines makes it possible to increase the system efficiency by 
35-49% with respect to the simple gas turbine cycle. 


COMBINED MARINE DIESEL ENGINE/ 
STEAM TURBINE SYSTEMS 


The calculations of the combined Diesel engine/steam 
turbine systems took into account two marine low-speed 
engines: 9RTA96C made by Warstsila [6] and 9K98MC made 
by MAN Diesel & Turbo [5] which revealed similar power 


outputs as the gas turbines. The analyses were performed in 
the combined system with the steam turbine cycle with a single 
pressure boiler, according to a scheme shown in Fig. 3, and the 
two-pressure boiler, see Fig. 4. The calculations were performed 
based on the producer’s data for the reference point according 
to ISO Conditions: ambient air temperature equal to 25°C and 
barometric pressure equal to 1 bar, Table 3. 


Tab. 3. Parameters of marine low-speed Diesel engines 
used in the combined Diesel engine/steam turbine system [ 5, 6 ] 


bep 


Type 
of engine 


kg/s I 


9K98MC 48762 | 134.25 | 232.8 | 0.482 | 174.9 
9RTA96C | 46332 | 104.50 | 271.0 | 0.506 | 166.8 


In that case the steam turbine system was analysed for the 
feed water temperature equal to tpw= 120°C. The results of the 
combined cycle calculations are given in Table 4. Like for the 
combined system with gas turbines, the caloric value of the 
fuel was assumed equal to Wu = 42700 kJ/kg. 

The power of the steam turbine in the combined system 
ranges from 5.3% to 6.8% of the piston engine power. The 
steam turbine power is larger for the steam turbine cycle with 
the two-pressure boiler, but this increase in not significant as 
compared to the power of the turbine with the single—pressure 
boiler. The turbine power increases with the increasing 
temperature of the piston engine exhaust gas. The pressure of 
the live steam for the two-pressure cycle is about 2.5 times as 
high as that for the single-pressure cycle and it increases with 
the increase of the temperature t, of the engine exhaust gas. 


Tab. 4. Combined Diesel engine — steam turbine cycle 


9K98MC 9RTA96C 


46332 


104.504 


2.1467 


271 


Single-pressure | Two-pressure 


Nsr 


Nsr/N D 


Neombi 


be .ompi 


m,/Np 


m,/m, 


m,,/m, 


m,/Np 
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The live steam mass flow rate for the single-pressure system 
is larger than that for the two-pressure system. When related to 
the engine power, this parameter differs insignificantly for the 
two examined engines. The extraction steam mass flow rate mg, 
as related to that of the live steam, is almost twice as small in 
the single-pressure system as in the two-pressure system. In the 
two-pressure system the mass flow rate m, of the low-pressure 
steam which feeds the steam turbine is equal to 59 + 73 % of 
the live steam mass flow rate and increases with the increase 
of the piston engine exhaust gas temperature. 

In the combined system the cycle efficiency increases by 
5.6 + 6.7% with respect to the power plant with only a Diesel 
engine. 


COMPARISON OF COMBINED GAS 
TURBINE/STEAM TURBINE AND MARINE 
LOW-SPEED DIESEL ENGINE/STEAM 
TURBINE SYSTEMS 


The use of combined systems with the gas turbine or 
the Diesel engine cooperating with the steam turbine cycle 
increases the power of the system, and simultaneously increases 
the power plant efficiency. 

In the systems with gas turbines the efficiency increase with 
respect to the simple gas turbine cycle (35 + 49%) is much 
larger than that observed in the combined cycle with the Diesel 
engine, the increase of an order of 5.5 + 6.7% with respect to 
the Diesel engine power, see Table 5. 


Tab. 5. Comparing combined cycles 


Combined gas turbine/steam turbine system 
CYCLON | 570 | 49.03 | -32.91 149.0 
GTX100 546 | 40.76 | -28.96 140.8 
GT10B 
LM2500 
GT35 


Combined Diesel engine/steam turbine system 


9K98MC 
9RTA96C 


232.8 
271 


5.58 
6.68 


-5.28 
-6.36 


105.6 131 
106.8 131 


The two combined systems reveal comparable overall 
efficiencies, and the resultant specific fuel consumption is also 
similar. It is noteworthy, however, that the heavy fuel burned 
in the marine low-speed engines is cheaper than the marine 
Diesel oil burned in gas turbines. This, in combination with 
comparable fuel consumption rates, results in lower fuel costs 
of the combined system with the Diesel engine. 

Table 6 gives the energy balance of the combined system 
with the gas turbines or the marine low-speed Diesel engine. 

In the energy balance of the combined system, the power 
of the main engine is equal to 35%, on average, for the gas 
turbine and 49% for the Diesel engine in relation to the energy 
delivered in the fuel to the engine, while the steam turbine 
power is equal to 15% for the system with the gas turbine and 
3% for the system with the Diesel engine. The stack loss and 
the condenser loss amount to 15% and 33% respectively in 


the system with the gas turbine, and to 15% and 11% in the 
systems with the Diesel engine. The comparable levels of the 
stack losses for these two systems result from their comparable 
exhaust gas temperatures. 

Other losses in the system with the Diesel engine amount to 
about 21%, as compared to 3% in the systems with gas turbines. 
In the Diesel engine the other losses include the heat taken over 
in the scavenge air coolers - about 15%, in the lubricating oil 
cooler — about 3%, and in the jacket water cooler 4%. 


Tab. 6. Energy balance of the combined cycle 


Type of engine 

Energy input 

Engine Power 
Steam Tur- 
bine Power 
Condenser 
Stack Loss 


SC 
xe 
xe 
xe 
xX 


CYCLON 
GTX100 
GT10B 


LM2500 
GT35 
9K98MC | 100 | 48.20 | 2.54 | 10.77 | 16.50 | 21.99 
9RTA96C | 100 | 50.55 | 3.17 | 11.42 | 15.68 | 19.18 
9K98MC | 100 | 48.20 | 2.69 | 11.38 | 15.65 | 22.07 
9RTA96C | 100 | 50.55 | 3.43 | 13.16 | 13.58 | 19.28 


Steam turbine cycles used in the combined system 
composed of gas turbines should include a two-pressure boiler, 
the least, to allow good utilisation of the energy of the high- 
temperature exhaust gas. 

In the combined systems comprising a low-speed Diesel 
engine, a single-pressure boiler can be used in the steam turbine 
cycle, as it can secure good utilisation of the exhaust gas heat 
(much lower temperatures of the exhaust gas, compared to 
those recorded in the gas turbine). 

The live steam parameters in the combined system with 
gas turbines are much higher than those in the system with the 
Diesel engine. 

The power outputs of the steam turbines in the systems 
with gas turbines are many times higher than those in the 
systems with the Diesel engine, compare Tables 2 and 4. For 
comparable powers of the combined system, the power of 
the steam turbine amounts to 26 + 33% of the power of the 
entire power plant for the systems with the gas turbines, and 
to 5.3 + 6.9% in the systems with the Diesel engine. In those 
latter systems the power output of the Diesel engine amounts 
to 94% of the power plant power, while in the systems with the 
gas turbines the power of the gas turbine amounts to 67 + 74%. 
For the same total powers of the combined system the power 
of the Diesel engine is larger 1.3+1.4 times than that of the 
gas turbine. 

The exit temperatures texy of the exhaust gas in the 
combined system for the cases with the gas turbine or the Diesel 
engine are comparable, Table 5. 
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FINAL CONCLUSIONS 


There is a possibility to use a combined system composed 
of the Diesel engine or a gas turbine as the leading engine, and 
the steam turbine cycle which utilises the heat contained in 
the engine exhaust gas. These systems reach thermodynamic 
efficiency comparable with the combined systems of gas 
turbines cooperating with steam turbines. 


Power outputs of the combined systems 

Depending on the adopted variant and main engine load, the 
use of the combined system can increase the power output of 
the power plant by 35 + 49% in the systems with gas turbines, 
and by 5 + 7% in systems with Diesel engines, with respect to 
a conventional power plant for the same fuel mass flow rate. 
Additional power output is generated due to recovery of the 
energy contained in the exhaust gas from the piston internal 
combustion engine or the gas turbine. This way the combined 
system reduces the specific fuel consumption by 26 + 36% 
for cycles with gas turbines and by 5 + 6.4% for the cycle 
with a Diesel engine, with respect to the conventional power 
plant. 


Efficiencies of the combined systems 

The use of the combined system for ship propulsion 
increases the system efficiency, which leads to the reduction 
of the specific fuel consumption and, additionally, increases 
the power of the of propulsion system, without additional 
fuel consumption. Like the power output, the efficiency of the 
combined system increases with respect to the conventional 
power plant, reaching the level of about 43 + 54% for maximal 
power outputs. The efficiency levels for the combined gas 
turbine/steam turbine and Diesel engine/steam turbine systems 
are comparable. 


Overall power plant dimensions 

The power plant with the combined system with gas 
turbines is smaller both in mass and dimensions than that with 
the Diesel engine. 


Economic aspects 

The technical and economic analysis of power plant 
operation should be performed to justify the use of such a power 
plant. The analysis should take into account both investment 
and operating costs. At comparable fuel consumption in both 
combined systems, the cost of burning the heavy fuel oil in 
the marine power plant with the Diesel engine is lower, as 
the combined systems with gas turbines burn more expensive 
marine Diesel oils. But in case of sailing in zones with limits 
imposed on heavy fuel burning, the use of a combined power 
plant with gas turbines can be more economically justified. 


Depending on the adopted solution, the use of the combined 
power plant makes it possible to reach the assumed power 
of the propulsion system at smaller load of the main internal 
combustion engine, and to reduce the fuel consumption. 
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NOMENCLATURE 

be — specific fuel consumption 

m — mass flow rate 

m,; | — mass flow rate of the low-pressure steam in the turbine 

N — power 

No — power of the gas turbine or of the piston internal 
combustion engine 

p — pressure 

Pi — low-pressure steam pressure 

ti — temperature of the exhaust gas from the piston internal 
combustion engine or the gas turbine 

n — efficiency 

SUBSCRIPTS 

combi — combined system 

D — deaerator 

D — piston internal combustion engine 

exh — outlet duct 

F — fuel 

FW — feed water 

GT — gas turbine 

K — parameters in the condenser 

o — live steam, calculation point 


ST — steam turbine 
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Exhaust gas temperature measurements 
in diagnostic examination of naval 
gas turbine engines 


Part III 
Diagnostic and operating tolerances 


Zbigniew Korczewski, Prof. 
Gdansk University of Technology 


ABSTRACT 


The third part of the article presents a method for detecting failures of the automatic 
engine control system with the aid of an exhaust gas temperature setter, specially designed 
and machined for this purpose. It also presents a procedure of identifying the operating 
tolerances and determining the diagnostic tolerances for the exhaust gas temperature 
recorded in the naval turbine engine during the start-up and acceleration processes. The 
diagnostic tolerances were determined using the statistic inference, based on the hypothesis 
about the normal distribution of the starting exhaust gas temperature dispersion at the 


initial time of engine operation. The above hypothesis was verified using the non-parametric statistic test 

X for examining the consistency of the empirical distribution with the assumed normal distribution. As 

a result of the examination, satisfactory convergence of the compared distributions was obtained which 

made the basis for assuming the three-sigma limits of the diagnostic tolerance for the analysed engine 
control parameter. 


Key words: technical diagnostics; naval turbine engines; exhaust gas temperature; 
operating and diagnostic tolerances 


INTRODUCTION 


The general assessment of the technical state of a naval 
turbine engine can be done based on the values of basic 
parameters, such as power, specific fuel consumption, rotational 
speeds of rotating units, and/or mass flow rates of the delivered 
fuel or the thermodynamic medium flowing through the 
engine. All these parameters characterise, in a general way, 
the quality of engine operation [1, 2, 3, 6]. However, in case 
of most serial engines, the above parameters cannot be directly 
measured (their measurement is extremely complicated and/or 
economically unjustified) in conditions of engine operation 
in the naval power plant. Moreover, these parameters do not 
deliver sufficiently precise information on the real course of 
the physical processes taking place in basic functional modules 
of the engine (flow part, kinematic system, fuel installation, 
automatic control system, etc.) [4, 5, 7, 8, 14]. That is why 
the operational diagnostics of turbine engines makes use of 
a set of auxiliary control parameters which indirectly define 
the basic parameters and, additionally, provide opportunities 
for indicating the areas in which the largest energy losses are 
observed, as well as for identifying and localising known 
and recognisable failures of engine sub-assemblies and 
constructional elements during engine operation. 


The real value of each control parameter differs from its 
calculated value, which results from differences in machining 
(inaccuracies of the realised technological process), the action 
of external agents (temperature, humidity, dust, vibrations, 
sea undulation, etc.), impurities, and/or wear and tear of 
constructional elements (reversible and irreversible changes of 
the surface layer) [7,8,9]. These actions are of random nature 
(random events) and, as a consequence, generate random 
changes of values of the structure parameters [2, 7]. Their 
effect on the technical state of the engine can also be observed 
as changes of the probability distribution of the occurrence of 
these parameters — the expected value and the average deviation 
(dispersion around the mean value for the set of engines in 
operation). The smaller the changes of the expected value and 
the mean deviation of the control parameter, the higher the 
correctness of engine operation and, consequently, its lifetime 
and reliability. 

The correctness of operation of a constructional module or 
unit, as well as that of the entire turbine engine, is defined by 
the tolerances of parameters which reflect the most desirable 
and permissible courses of the realised physical processes: 

e Operating tolerances (maximal permissible) — determined by 
the producer (designer) based on the tests done on prototype 
copies. They define precisely the maximum ranges within 
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which the values of the control parameters are to be kept 
during the engine operation process to secure that the engine 
performance parameters do not drop below permissible 
levels, the engine operation is undisturbed and reliable, 
and that the technical durability is good and in line with 
the limits of the operating potential for particular hours and 
days. When the value of any control parameter exceeds the 
set operating tolerance limits, it is a signal of unacceptable 
disturbance of energy conversion processes realised in the 
engine, which threats with its failure (unstable operation of 
the compressor, for instance) [3, 7, 8, 11, 12, 14]. 

e Diagnostic tolerances — determined via statistical inference 
of the results of the tests done on a sufficiently numerous 
set of new and correctly adjusted engines (revealing full 
technical ability), compared with the results of operating 
tests of the same engines is successive stages of their use, 
provided that they still maintain a good functioning level, 
which is assessed from changes of basic parameters. During 
the time of engine operation the values of the control 
parameters shift towards the diagnostic tolerance limits, 
which signals the symptoms of small changes of engine’s 
technical state, characteristic for the “approaching” state of 
inability. At the same time, the above changes of the basic 
parameters reflect the dependence of the technical state of 
the engine on the time of its operation, which manifests 
itself as changes of the expected values of the control 
parameters and their average deviations. 


Consequently, a conclusion can be formulated that the 
tolerances of the control parameters for diagnostic purposes 
should be much smaller (more restrictive) that the operating 
tolerances. 


IDENTIFYING FAILURES OF THE 
AUTOMATICS SYSTEM (OPERATING 
TOLERANCES) 


Some most typical operating failures of the engine 
automatics system can be identified using simulation analyses 
done with the aid of specially designed setters (testers) on 
the non-operating engine [15]. These devices control the 
correctness of operation of the engine protection system against 
the excessive increase of the exhaust gas temperature. The tests 
aim at checking correct threshold settings for the (maximal) 
exhaust gas temperatures behind the exhaust gas generator, and 
the temperature grow rates at which the final control units in 
the automatics system open: 

- the fuel overflow valve (ZPP), 
- the main fuel valve (ZGP), 
- the immediate (emergency) engine switch-off valve 

(ZA). 


The values determined in the above way are considered the 
limits of the operating tolerance zone for permissible changes 
of parameters in the process of engine start-up and acceleration. 
Sample results of these tests are shown in Fig. | in the form of 
simulated dynamic time-histories of exhaust gas temperature 
changes. To assess the ability of the examined engine, the time- 
histories obtained from the simulation test are compared with 
the real curves recorded during earlier start-ups (see: part II 
— Unsteady states). 

The comparison analysis of the results recorded during 
many years of operation of the UGT type ZORYA naval engines 
which were the object of diagnostic supervision [15] has made 
it possible to recognise characteristic thresholds for initiating 
certain actions of the protection system. These thresholds are 
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Fig. 1. Simulated (dashed lines) and real (solid lines) time-histories of the 
exhaust gas temperature recorded behind the low pressure turbine during 
diagnostic start-up tests of the UGT type ZORYA engine. 
1- input function dT}, /dt = 112 K/s - ZPP(-), ZGP(+), ZA(+); 
2 — input function dT /dt = 88 K/s - ZPP(-), ZGP(+), ZA(+); 
3 — input function dT}, /de = 84 K/s - ZPP(+), ZGP(+), ZA(+); 
4 — input function dT}, /de = 48 K/s - ZPP(-), ZGP(-), ZA(-); 
5 — successful engine start-up aT,» /dt = 64 K/s ; 
6 — unsuccessful start-up dT, /dt = 82 K/s; 


A — automatic engine switch-off 


considered the limits of the operating tolerance zone for the 

time derivative ofthe averaged exhaust gas temperature behind 

the exhaust gas generator: 

° dTi jdr > 80 K/s- the engine is automatically switched-off 
- ZPP(+), ZGP(+), ZA(+); 

e dT,,/dt = 50+ 80 K/s — the engine can be started with 
automatic fuel dosage - ZPP(+), ZGP(-), ZA(-); 

e dT,,/dt< 50 K/s — the engine is starter automatically 
without the action of the protection system - ZPP(-), ZGP(-), 
ZA(-). 


Having analysed the data shown in Fig. 1 we can conclude 
that the system of thermal engine protection is out of order 
in this case. The fuel overflow valve (ZPP) does not work, 
irrespective of the rate of exhaust gas temperature grow (curves 
1 and 2). After a detailed analysis of the recorded symptoms of 
engine inability, the damaged engine automatics elements were 
localised and replaced. The next simulation tests confirmed the 
correctness of the formulated diagnosis, which is illustrated in 
Fig. | by the line 3 recorded after removing the defect. 

It also results from the performed tests that the system 
has the zone of insensitivity ranging between 273 and 593 K, 
irrespective of the dynamics of exhaust gas temperature 
grow. 


METHOD OF DEFINING DIAGNOSTIC 
TOLERANCES 


Measurements of the control parameters which characterise 
the quality of the turbine engine operation in the processes 
of start-up, acceleration and deceleration of rotor units, the 
operation at steady-load ranges, and the engine switch-off are 
done in the very initial stage of engine use, directly after its 
installation in the power plant room [14, 15]. Discrete values 


of the diagnostic parameters recorded in these measurements 
are modelled as continuous and one-dimensional random 
variables with the certain probability distribution, expected 
value, and variance being the measure of the dispersion of 
the measured results. In this situation the diagnostic tolerance 
calculation method is mainly limited to finding the form of the 
distribution of the dispersion of the observed discrete values of 
the diagnostic parameter around the average value for the entire 
population of engines in operation, and then to determining the 
parameters of this distribution. One of methods making use 
of the statistic inference to find the form of the distribution 
of the random variable consists in verifying the hypothesis 
that the examined distribution of the discrete values of the 
diagnostic parameter can be approximated using an assumed 
theoretical distribution which is described by a known random 
variable probability density function at time t = 0 [2, 10]. Once 
the form of the function of the theoretical distribution and 
its characteristic parameters (the expected value m and the 
average deviation o) are known, we can determine the interval 
of tolerance for the values of the diagnostic parameters, for 
instance m + 30 for the normal distribution. 

A method for determining the limits of the diagnostic 
tolerance zone for the rate of grow of the exhaust gas 
temperature behind the exhaust gas generator in the start-up 
process), dT,,/dt, which take into account the unrepeatability 
of production and the effect of external operating agents, will 
be demonstrated on the population of 46 UGT3000 type Zorya 
turbine engines which were introduced to operation in naval 
power plants [14,15]. Table 1 collects a set of realisations of 
the numerical values of the starting temperatures recorded in 
the start-up process of the new engines (or directly after general 
overhaul) in the state of full technical ability (the values are 
reduced to normal atmospheric conditions) [3, 13]. 


Tab, 1. Set of realisations of numerical values of the starting temperature 
dT, /dt (random variable X(t)) recorded during the start-up in the initial 
time of operation of UGT3000 engines installed in naval power plants 


Assuming that the distribution of the values of the analysed 
diagnostic parameter in time t = 0 is normal, which is confirmed 
by diagnostic tests of aircraft engines [2], a zero hypothesis 
H, is formulated that the distribution function of the empirical 
distribution shown in Table | is consistent with the normal 
distribution function. The zero hypothesis was verified using 
the non-parametric statistical test y2, which does not require 
the information on normal distribution parameters, i.e. the 
expected value m and the variance o? [10]. In this situation the 
statistical estimation of the expected value and the variance for 
the calculated statistics is made based on the arithmetic mean X 
and the standard deviation © of the measured results. 

For testing purposes, the interval within which the values 
of the random variable are included X(t) = {x;} was divided 
into six equal sub-intervals of 3,0 K/s in length (the widths 
of the intervals do not have to be equal) — Table 2. The limits 
of the intervals were chosen in such a way that the numbers 
of the results of measurement in successive intervals were 
not excessively small (5-8 results in each interval as the 
minimum). 


Tab. 2. Scheme of y} statistics calculations 


Number of results 
of measurement 
in i-th interval 


Limits of 


Number of ; 
intervals 


interval 


Theoretical probability 


Normalised 
limits of 
intervals 


Value of Laplace function 
for beginnings of 
intervals 
®p(zi) 


0.1065 
0.3507 
0.4687 


Value of statistics for i-th 


to obtain the result 
from i-th interval 


Pi = ®, @.1)- ®, (;) 


Number of 
interval 


Theoretical number of 
results in i-th interval 


n’ pi 


interval 


TE 2 (n; =n-pi) 
n-p; 


0.6668 
1.2349 


6.8678 0.1097 


46 2.8012 
e arithmetic mean X = 18.97826, « standard deviation © = 3.784919 


D) Further in the text, more briefly referred to as the “starting temperature”. 
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If the number of the results of measurement in the interval 
is smaller than 5, this interval is to be linked with one of the 
neighbouring intervals [10]. 

In successive steps we calculate: 

e the normalised limits of the intervals with respect to X, re- 
calculated to © units: 


z= (1) 


The beginning of the first interval is -00, while the end of 
the last interval is +0. The end of the previous interval is 
simultaneously the beginning of the next interval; 

e the normalised Laplace function of the normal distribution 
for the beginnings ®p(z1) of the intervals — from statistical 
tables [10]; 

e the theoretical probability to obtain the result from the i-th 
interval: 


P:i = PZ) — ®,(Z;) (2) 
e the theoretical number of results in the i-th interval: n - p;; 
e the y2 statistics for each interval: 


2 
i= (nj —n-pj) (3) 
l n pi T 
The last step includes calculating the value of the statistics 
for the analysed set of random variable realisations: 


x? =}_x; =2.8012 (4) 
i=l 


The calculated value of the y,? statistics should be compared 
with its critical value determined from the statistical tables [10], 
at the assumed significance level a and the number of freedom 
degrees f calculated from the following relation: 


f=l1-k-1 (5) 


where: 
| = 
k= 


number of intervals in the distributive series, 
number of estimated parameters of the verified 
distribution. 


Two parameters characterising the normal distribution: 
X and 6, were evaluated based on the empirical data. For the 
number of degrees of freedom: f = 2 (1 = 5, k = 2) and the 
assumed significance level: a = 0.05 the critical value of the 
statistics which was read in the %? distribution tables is equal 
to %2 = 5.991 [10]. 

The calculated value of the statistics is equal to y,2= 2.8012 
and is smaller than the critical value y,,2(a = 0.05, f=2)=5.991. 
That means that there is no ground for rejecting the zero 
hypothesis and, consequently, that the analysed empirical 
distribution can be considered consistent with the normal 
distribution. Therefore the probability density function of the 
random variable X(t) = {x,} in the initial operation stage (time 
t= 0) can be given by the formula: 


f(x) = - exp|— (6) 


1 
o-V2n 

Having known the formula for the distribution and its 
basic parameters we can calculate the limits of the diagnostic 
tolerance zone assuming the value interval m+3o as the 
diagnostic limiting conditions for the control parameter of 
interest. For the analysed starting temperature the calculated 
expected value was equal to m = X= 18.98 K/s and the average 
deviation was o = © = 3.78 K/s, which gives the following 
diagnostic tolerances of the analysed control parameter: 
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7.64 < x < 30.32 (7) 


It results from the performed calculations that the diagnostic 
tolerances of the rate of exhaust gas temperature grow in the 
start-up process of the turbine engine of UGT3000 type are 
limited from one side, which means that the limiting value of 
this parameter assessed for diagnostic purposes must not exceed 
30 K/s in a new engine. 

Taking advantage, in turn, of the symmetry of the normal 
distribution and the values of the Laplace function ®,(z;) we 
can calculate the probability that the random variable X(t) takes 
a value from within this interval: 


P{X e (m +30)} = ®, (3)-@, (-3) = 0.9973 (8) 


Thus, assuming the three-sigma limits for the diagnostic 
tolerances of the analysed control parameter means that less 
than 3 ə (3 out of 1000 as the maximum) of possible results of 
measurement differ from the expected value by more than 30. 

In the process of engine operation on a ship, the action 
of unfavourable operating agents leads to ageing, wear, and 
pollution of constructional elements. These processes are 
permanent and irreversible, and they always accompany the 
operation of internal combustion engines in marine conditions. 
All this leads to the decrease of the usable potential, the measure 
of which are the deformations of the shape of probability 
density function defined by the changes of the expected value 
and the average deviation, while the normal distribution of the 
random variable (dispersion of the analysed control parameter) 
is still being preserved. As a consequence, the dispersion of the 
exhaust gas temperature grow rate in the engine start-up process 
increases, and the limits of the diagnostic tolerances, which 
under no circumstances can exceed the operating tolerances, 
become wider — Fig. 2. We accept that with the increasing total 
time of engine operation its functioning is getting less correct, 
but still preserves all requirements of safe use. 
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Fig. 2. Deformations of the shape of the probability density function 
for the analysed random variable vs. the total time 
of operation of the turbine engine 


It results from the numerical data shown in Fig. 1 that after 
calculating the shapes ofthe probability density function ofthe 
random variable in successive stages of engine operation (for 
instance, in annual control and measurement cycles, or after 
each 500 hours of operation) we can determine the trend lines 
for the changes of the expected value and the upper tolerance 
limit of the set or realisations of the diagnostic parameter. We 
can also assess the horizon for the prognosis of correct engine 
operation (without failures) until time t,,, based on the crossing 
point of the line of upper diagnostic tolerance limit with the 
line of the operating tolerance limit. For this purpose the value 


of the diagnostic parameter is to be measured on the examined 
engine - for instance, the point: 


* 
t=tegl 7 
t=tegl 
and the upper diagnostic tolerance limit is to be extrapolated 
until it crosses the line of the operating tolerance limit 
(time t.) etc. This way we formulate the strategy of engine 
operation based on its current technical state, the reconstruction 
of which (adjustment, repair, replacement of an element) 
depends on individual features of constructional structure 
degradation, mainly connected with the nature of changes of 
engine load (number of start-ups, stop-down, accelerations 
and decelerations of rotor units, manoeuvres of the separate 
power turbine, etc.). 


CONCLUSIONS 


The presented sample application of the statistical method 
for determining the diagnostic tolerances of the rate of grow of 
the exhaust gas temperature recorded in the start-up process of 
a naval turbine engine can be extended onto remaining control 
parameters recorded during engine operation. A necessary 
condition for statistical inference is to adopt the assumption that 
in the initial stage of operation the examined engines maintain 
correct functioning, which is assessed from the values of basic 
parameters. In this situation the value of the observed control 
parameter shifts towards the limit of the diagnostic tolerance 
zone, thus generating a symptom of small changes of the 
technical state, characteristic for the “approaching” inability 
state. Simultaneously it reflects the dependence of the technical 
state of the engine on the time of its use, which manifests itself 
by changes of the expected values of the control parameter and 
its average deviations in successive engine operation stages. 
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Influence of pitting corrosion on strength 
of steel ships and offshore structures 
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ABSTRACT 


The present paper deals with the influence of pitting corrosion on mechanical properties 
of mild and low alloy steels and strength of steel structures under static and quasi-static 
loads. Pitting corrosion is avery important phenomenon that influence local strength 
of ship hull members. The present paper is based mainly on Japanese publications. 
Analyzing the standard tensile diagrams one can see that that pitting corrosion reduces 
the load corresponding to yield stress (YS) and, even more markedly, tensile strength UTS, 
reduces almost to zero plastic flow strains at YS load level and dramatically reduces the 


total elongation to fracture. Nominal UTS load for uniformly corroded specimens is higher than that for 
pitted specimens at the same average thickness loss. The strength related to the true fracture surface area 
(i.e. reduced by pits, taking into account real path of the crack through the pits) is almost independent of the 
thickness loss and is almost the same as for uniform thickness loss. For wide specimens the tensile strength 
depends mainly on the local deformation ability, and the maximum loading ability for large members, 
predicted on the base of the total true fracture surface area, can be overestimated. Concept of equivalent 
thickness loss for plates under bending and compression and for more complex structural models as well 
as relation between the average and the equivalent thickness losses has been presented. Approach based 
on degree of pitting (ratio of pitted area to total plate surface area) has been shown as a real and more 
convenient than the approach based on the thickness loss. 


Keywords: pitting corrosion, strength, shipbuilding steels, ship structures 


INTRODUCTION 


There are two groups of actions that can lead to damages 
of real structures: 
e Chemical or electrochemical, if undesirable, are identified 
with corrosion 
e mechanical, usually identified with stresses or strains 


Interaction of corrosive and mechanical factors is the most 
general case. Almost all corroding structures are stressed. 
Almost all stressed structures are exploited in any environment 
that is not neutral in mechanical damage process. Pure 
mechanical or pure corrosive damage can be considered as 
specific and unusual cases. In case of vacuum or eventually 
inert gas environment we can assume that damage is purely 
mechanical in nature, while for very low stress levels we can 
assume that damage is purely corrosive. 

The above is especially true for ships and offshore structures 
that are the main object of the present paper. Corrosive 
environment of these structures is also the main source of 
service loads, i.e. static or dynamic pressure of sea water. 
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The present paper concerns the influence of pitting 
corrosion on strength of steels and steel structures under static 
or quasi-static loads. The mild and low alloy steels used for 
welded marine structures are considered as insensitive to stress 
corrosion cracking in marine environment. Therefore only the 
pit nucleation and growth process should be considered as 
a complex mechanical-electrochemical phenomenon. Fracture 
process occurs in a relatively short time and influence of 
corrosive environment on this process can be considered as 
negligible. Pitting corrosion incubation and growth and the 
role of strain and/or stress in this process has been discussed 
elsewhere [1, 2]. The present paper is focused on the influence 
of pre-existing pits on mechanical properties of mild and low 
alloy steels and strength of the steel structures. 

Previous pitting corrosion influences the fatigue strength 
and life by the following ways: 

I reduction of the average material thickness that influence 
the nominal stress in the structure; 

II pitting-corrosion-induced roughness of the material surface 
that leads to local stress concentrations at some locations 
at the surface. 


The role of pitting in fatigue of materials and structures is 
commonly appreciated. The problem has been discussed in [2] 
on the wide base of references concerning different materials. 
In common opinion, pits, if they are present at the material, are 
almost always the potential sites of the fatigue crack initiation, 
thus they reduce the fatigue life. If strength of steel structures 
under static or quasi-static loadings is considered, the opinion 
is not so unanimous. 

Melchers [3] has formulated opinion that pitting corrosion 
has very little influence on the strength of well maintained 
structures, such as ships — general corrosion (loss of thickness) 
is more important. This opinion has been accepted by Wang 
et al. [4]. These authors have elaborated models for the time- 
related growth of the maximum depth of pitting because, in 
their opinion, the main menace is perforation of the structure 
walls and it can be dangerous especially for pipelines and tanks. 
However, extensive studies of the influence of pitting corrosion 
on the behaviors of materials and structures under static 
loading has been carried out in Japan by Nakai, Matsushita, 
Yamamoto and, in some papers, Arai [5 - 10], and they have 
shown that the influence of pitting on the steel strength is not 
negligible. General strength of ship hull is mainly affected by 
general corrosion that reduce thickness of shell and stiffeners 
therefore reduce midship section modulus. Pitting corrosion 
is a very important phenomenon that influence local strength 
of ship hull members. 

International Association of Classification Societies (ACS) 
[11, 12] and Polish Register of Shipping (PRS) [11, 12] have 
defined the required minimum thickness of plating in pits or 
other local corrosion areas (is to be greater than 70-75% of 
the as-built thickness), or average thickness across any cross 
section in the plating (is to be greater than the renewal thickness 
for general corrosion). These requirements are surely based on 
unknown considerations of strength of pitted steel. 


TENSILE DIAGRAM FOR PITTED STEEL 


The mentioned Japanese authors tested only conical pits 
that are observed on hold frames of bulk carriers. The pits were 
generated artificially by spot drilling of the both steel surfaces. 
Different pit patterns were manufactured with different values 
of degree of pit intensity DOP (ratio of pitted area to entire area 
of the specimen). In the Rules of Classification Societies DOP 
is called simply the pitting intensity. In order to verification 
of validity of such artificial pits testing, the results were 
compared to the results for actually corroded specimens [9]. It 
has appeared that specimens with orderly located artificial pits 
which had a fixed diameter equal to the average pit diameter in 
actual corroded specimens could simulate the actual corroded 
specimen where the pit size and distribution are random. Results 
of standard static tensile tests of pitted and smooth specimens 
are described in [6, 7, 9]. Exemplary tensile diagrams (load 
versus displacement) for small specimens are shown in Fig. 1 
[7]. It is evident that pitting: 

(i) reduced the load corresponding to yield stress; 

(ii) reduced almost to zero plastic flaw strain at the level of 

yield strength; 

reduced load, corresponding to ultimate tensile strength, 

even stronger than yield stress; 

markedly reduced total elongation to fracture; 

(v) total elongation for 13mm and 11mm uniformly thick 
specimens is the same. 


The same investigations [7] have shown that the effect 
of pitting depends on the specimen thickness, i.e. the 
nominal tensile strength and total elongation dropped more 
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Fig. 1. Results of small tensile specimens testing: a) view of specimens after 
testing; b) tensile diagram “load — displacement” [7] 
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significantly for higher pit-depth-to-plate-thickness ratios. 
For thinnest plated (10 mm) total elongation dropped from 
30% to about 5%. 

Matsushita et al [6] realized tensile tests of specimens with 
different DOP from a single one-side pit up to more than 80% 
of two-side surfaces covered by pitting. Pitting corrosion causes 
not only unevenness of the material surface but reduction of 
the average thickness too. The larger DOP leads to the larger 
average thickness loss. Results are plotted as a function of the 
thickness loss in Fig. 2 and Fig. 3. 

Trend lines for ideally uniform corrosion are also plotted 
in Fig. 2. All over the range of thickness loss values, the 
nominal strength for pitted specimens is lower than for 
“uniformly corroded” ones. The specimens fractured at the 
smallest area where the local thickness reduction was larger 
than the average value. Meanwhile, strength of the “uniform 
corrosion” specimen with a value of uniform thickness loss is 
always compared with the pitted specimen strength of the same 
average (not local) thickness loss. True fracture surface does 
not develop close to one surface approximately perpendicular 
to load axis, as that in non-pitted specimens, but jumps from 
one plane to another according to pits distribution. True fracture 
surface area can be evaluated as the area of the fracture surface 
before loading projected to the plane perpendicular to the load 
axis and practically can be substituted by the minimum cross- 
sectional area (i.e. locally reduced by pits) before loading 
perpendicular to the load axis [9]. In Matsushita et al opinion, 
the true strength that is based on the true fracture surface area 
is the same for uniform and non-uniform thickness loss (Fig. 
2b). In our opinion, however, the above explanation is only 
partially valid. If the only specimen of another type (wide 
specimen — type B) is not taken into account, the strength of 
pitted specimens for average thickness loss below 2 mm will 
be 527-600 MPa, while for average thickness loss above 2 mm 
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Fig. 2. Tensile strength of pitted steel versus thickness 
loss for type A (small specimens) and type B (wide specimens): 
a) nominal strength P/A, (P is the load, A, is the initial cross-section area); 


b) true strength P/A (A is the true area of fracture surface, 
i.e. reduced by pits) [6] 


the strength will be 400-527 MPa thus it is not independent of 
the intensity of pitting. In spite of that, the maximum load that 
can be carried by a specimen can be approximately predicted by 
the above approach based on true fracture surface or minimum 
cross-sectional area [6, 9]. 

The last conclusion concern so called small specimens. 
The situation is more complex in case pf large (much more 
wide) specimens. Fracture of small specimen occurs when 
plastic deformation develops throughout the minimum cross 
section. In much more wide specimens, having areas with pits 
of different size and distribution, the deformation ability differ 
from location to location. As a result, fracture occur locally at 
this regions where elongation reaches critical value first, i.e. 
where the elongation ability is the lowest due to local pits sizes 
and distribution, this local fracture is the dominant factor that 
determines the maximum load ability of the specimen modeling 
a structural member [9]. That is why the maximum loading 
ability for large members, predicted on the base of total true 
fracture surface area, is overestimated, while for small members 
this approach was successful [10]. 

Analogous diagram for total elongation is shown in 
Fig. 3. The uniform thickness loss does not cause reduction 
of total elongation (see also Fig. 1). Thus the observed very 
significant drop of values of total elongation has to be entirely 
attributed to non-uniform character of pitting corrosion. The 
drop is very rapid for small values of thickness loss (also 
values DOP are small). Then relationship saturates and the 
total elongation does not depend on thickness loss and DOP. 
Nakai et al [7] noted the minimum, and next even an increase 
of total elongation as DOP is increased further. This is perhaps 
due to decrease of the ratio of pit depth to average thickness 
loss. These data are for 10 mm thick material. For thicker 
materials the reduction of elongation is also evident but not 
so dramatic. 
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Fig. 3. Relationship between thickness loss 
and total elongation in tensile test [6] 


ULTIMATE STRENGTH 
OF STRUCTURAL MODELS 


Lateral-torsional buckling of whole hold frames, and local 
buckling and fracture of web and face plates in the lower parts 
of hold frames, is damage typically seen in the hold frames of 
bulk carriers [5]. Therefore extensive program of FE analyses 
has been realized [5] with square plates containing conical 
corrosion pits (aspect ratio depth-to-width of pit equal to 
0.125) of different size and distribution. Compressive load plus 
bending moment have been applied in different proportions: 
from pure compression (stress gradient equal 0) to pure bending 
(gradient equal 2). 

As for tensile test, for these loads it also has appeared that 
ultimate strength of pitted plates is smaller than or equal to that 
of uniformly corroded plate for the same average thickness 
loss. There is a tendency for the strength reduction to increase 
as the pit diameter and depth is increased. Thus a prediction 
of the strength of pitted plates using the average thickness loss 
would lead to non-conservative results. However, the average 
thickness loss or maximum pit depths are not always dominant 
factors for ultimate strength but the strength depends also on 
type of pit distribution on the both sides of the plate [5, 9]. 

Nakai and co-workers analyse equivalent thickness of plate 
for ultimate tensile strength (t,) that is defined as the thickness of 
an uniformly corroded plate whose ultimate strength is the same 
as that of pitted plate. Value of the equivalent thickness is almost 
independent of the loading mode: from pure compression, 
through different combinations of compression and bending, 
up to pure bending [5]. Nakai and co-workers have compared 
these values of equivalent thickness to the values evaluated by 
tensile test, where the following proportion is valid: 


te/to = O, / Ouo (1) 


where t, is the original plate thickness, o, and O, are the 
nominal UTS for the plate with and without pits, respectively. 
The comparison has been made for the same degree of pit 
intensity (DOP). It has appeared that equivalent thickness 
evaluated in tensile test of pitted specimens is smaller than that 
for evaluated for plates subjected to in-plane compression and 
bending, especially for the thinnest specimens. Thus tensile test 
represent the most dangerous case. If an evaluation method for 
the effect of pitting corrosion on the equivalent thickness (t,) 
under tensile loading is established, a reduction of t, of pitted 
plates under in-plane compression, bending and shear will be 
predicted conservatively by the same method. 

Equivalent thickness of plates with pitting corrosion (t,) for 
the tensile strength is almost equal to average thickness at the 
minimum cross section (t [10]: 


wenn) 


t, = retin (2 ) 


Equivalent thickness loss (t„) is defined as thickness loss of 
uniformly corroded plate that has the same tensile strength as 
a plate with pitting corrosion. It can be calculated by: 


ta = to z t. (3) 
Average thickness loss (t,,) of pitted plates is defined by: 
ta = to _ lie (4) 


where: 
t = 


ave 


the average thickness of the pitted plates. 


Finite Element (FE) analyses of pitted square plates 
under uniaxial compression, biaxial compression, shear and 
combinations if the above mentioned loadings showed that 
the equivalent thickness can conservatively be predicted by 
equation (3) and the following equation [10]: 


ty = 1.25 - ty (5) 


Thus ultimate strength of pitted plate of an average 
thickness (tve ) is lower than the strength of uniformly corroded 
plate of the same thickness. This means that unevenness of the 
plate surface due to pitting corrosion reduce ultimate strength 
of the plate. 

FE analyses of structural model: beam (web plate, face 
plate and shell) under 4-points bending has been done [10]. 
Web plates were pitted. It has appeared that for the following 
cases: (i) the ultimate strength is reached accompanying local 
buckling; (ii) the structural models collapse without buckling; 
(iii) shear buckling of web plates occurs; and (iv) web crippling 
occurs, equivalent thickness of the web plates can be predicted 
by equation (3) and the following equation: 


ta = 1.44: ta (6) 


The above equations mean that the influence of the pitting- 
corrosion-induced unevenness of the plate surface on the 
ultimate strength of more complex structural models is more 
detrimental than for simple square plates. 


APPROACH BASED ON DEGREE OF 
PITTING (DOP) oss THICKNESS 
L 


Nakai and co-workers conclude that for non-uniform 
random pitting, the average thickness loss of a whole plate, or 
for the minimum cross section, is often difficult to determine. 
Each pit is covered with a heavy rust blister which is difficult 
to remove, even with hammers. Even if the rust blisters and 
rust from non-pitted areas are removed, correct thickness 
measurement is difficult due to great unevenness of the surface. 
[9, 10]. The average thickness loss or the average thickness loss 
at the minimum cross section cannot be determined directly by 
observation of the plate surface. Meanwhile, degree of pitting 
(DOP) can be approximately evaluated by observation of the 
plate surface also without removing the heavy rust blisters, 
because Nakai and co-workers has shown that the diameter of 
the rust blisters and that of the corresponding pits are almost 
the same [10] and degree of blisters intensity (ratio of the 
surface area covered with rust blisters to the entire surface 
area) corresponds to DOP. Therefore, to elaborate a method 
of evaluating the tensile strength of pitted members using the 
value of DOP, is a task of primary importance. 

First attempt to solve this problem has been made by the 
same authors [10]. Diagrams of thickness loss of plates versus 
values of DOP are shown in Fig. 4. The diagrams include 


average thickness loss t,, and effective thickness loss t, for 
ultimate strength calculated for different cases: square plates 
under combined loading of compression and shear [equation 
(5)], more complex structural models (equation (6)) and for 
ultimate tensile strength. From ultimate-strength point of 
view, average thickness loss is not so important as effective 
thickness loss. Equivalent thickness loss for all three cases 
is increases almost linearly with the increase of DOP when 
DOP is smaller than approximately 75% [10]. The trend is 
especially evident foe equivalent thickness loss for ultimate 
tensile strength. Equivalent thickness loss for tensile strength 
(that equals to average thickness loss at minimum cross 
section) is larger than that for other loading conditions; only 
for largest values of DOP it is equal to 1.44t,, that is specified 
for structural models The same observation has been made for 
plates under in-plane compression and bending [5] described 
above. Therefore, it is possible to ensure the structural integrity 
if the average thickness loss at the minimum cross section 
is within the allowable corrosion level specified in rules for 
uniformly corroded members [10]. Allowable value of DOP 
could be evaluated on the base of allowable thickness loss 
for uniform corrosion, equivalent thickness loss for the pitted 
plate, and Fig. 4. 
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Fig. 4. Diagram of average thickness loss t„ equivalent thickness loss 
under combination of compression and shear (1.25t.,,), equivalent thickness 
loss for ultimate strength of structural models (1.44t,,) and equivalent 
thickness loss for tensile strength at minimum cross section (ta, min) 
versus degree of pitting (DOP) [10] 


Analyzing progress of pitting corrosion as a function of 
DOP, it has appeared [10] that average thickness loss and 
standard deviation of the thickness loss vary with small 
scatter band, while large scatter band is observed for extreme 
values like maximum pit depth or average thickness loss at the 
minimum cross section. When average thickness loss exceeds 
about 2 mm for one side of plate, then DOP equals about 100%. 
After this stage the form of corrosion tends to a change from 
pitting corrosion to general corrosion. 

The same authors [8] have investigated experimentally 
and analytically (FEM) strength of 4 point bent welded frames 
contained of web plate, face plate and shell (tensile load act 
on the face plate). Artificial pits were made over a separated 
areas of the web plate. Lowest strength had a frame with a band 
of pits close to the compressed shell, higher strength had the 
frame with a band of pits situated at the center of the web 
height, highest strength had the frame with pits close to tensed 
face plate. Pits situated at the center of the frame length (pure 
bending segment) more significantly reduced strength than 
that situated close to ends (high shear stress). It means that 
localization of pitting corrosion at the web plate surface is an 
additional factor influencing the strength of frame. Influence 
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of this position of pitting corrosion on the frame strength 
for longer, slender frames is less pronounced than for short 
frames. 


CONCLUSIONS 


1. Uniform thickness loss dos not influence the total elongation 
of specimens in tensile test. Pitting corrosion, even of 
small degree of pit intensity, very significantly reduce total 
elongation of the specimens. Pitting corrosion reduces the 
ultimate tensile strength of steels stronger than the uniform 
corrosion 

2. For tension and any combination of in-plane compression 
and bending, the ultimate strength of pitted plate is 
smaller than that of uniformly corroded plate with the 
same average thickness loss. There is a tendency for the 
strength reduction to increase with the as the pit diameter 
and depth is increased. Thus a prediction of the strength of 
pitted plates using the average thickness loss would lead to 
non-conservative results 

3. The equivalent thickness (t,), defined as the thickness of 
uniformly corroded plate whose ultimate strength is the 
same as that of pitted plate, for tensile loading is smaller 
than for any combination of in-plane compression and 
bending. Therefore it can be applied for conservative 
prediction of strength of members subjected to each of the 
loading modes. 

4. Itis possible and convenient to analyze the ultimate strength 
of pitted structures as a function of degree of pitting (DOP 
— the ratio of the total pitted area to the entire area of the 
plate) instead as a function of the plate thickness loss. 

5. Ultimate buckling strength of pitted plates and frames 
depends also of the position of the pitted areas and the type 
of pit distribution in every area. 
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The Ship Handling Research and Training Centre at Ilawa is owned by the Foundation for Safety of Navigation 
and Environment Protection, which is a joint venture between the Gdynia Maritime University, the Gdansk University 
of Technology and the City of Ilawa. 


Two main fields of activity of the Foundation are: 


Ə Training on ship handling. Since 1980 more than 2500 ship masters and pilots from 35 countries were trained at 
Itawa Centre. The Foundation for Safety of Navigation and Environment Protection, being non-profit organisation 
is reinvesting all spare funds in new facilities and each year to the existing facilities new models and new training 
areas were added. Existing training models each year are also modernised, that's why at present the Centre represents 
a modern facility perfectly capable to perform training on ship handling of shipmasters, pilots and tug masters. 


Research on ship's manoeuvrability. Many experimental and theoretical research programmes covering different 
problems of manoeuvrability (including human effect, harbour and waterway design) are successfully realised at 
the Centre. 


The Foundation possesses ISO 9001 quality certificate. 
Why training on ship handling? 


The safe handling of ships depends on many factors - on ship's manoeuvring characteristics, human factor (operator 
experience and skill, his behaviour in stressed situation, etc.), actual environmental conditions, and degree of water 
area restriction. 


Results of analysis of CRG (collisions, rammings and groundings) casualties show that in one third of all the 
human error is involved, and the same amount of CRG casualties is attributed to the poor controllability of ships. 
Training on ship handling is largely recommended by IMO as one of the most effective method for improving the 
safety at sea. The goal of the above training is to gain theoretical and practical knowledge on ship handling in a wide 
number of different situations met in practice at sea. 


For further information please contact: 
The Foundation for Safety of Navigation and Environment Protection 


Head office: Ship Handling Centre: 
36, Chrzanowskiego Street 14-200 ILAWA-KAMIONKA, POLAND 
80-278 GDANSK, POLAND tel./fax: +48 (0) 89 648 74 90 
tel./fax: +48 (0) 58 341 59 19 e-mail: office@ilawashiphandling.com.pl 

e-mail: office@portilawa.com 
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